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SUMMARY 


The  project  objective  was  the  development  of  fullwave  analyses  for  the  character¬ 
ization  of  microstrip  antennas  and  arrays  on  multi-layer  anisotropic  substrates. 
Many  other  researchers  have  addressed  the  problem  of  strip  antennas  on  isotropic 
substrates.  The  material  constants  are  very  important  to  the  performance  of 
printed  antennas.  In  practice,  many  dielectric  and  magnetic  media  are  anisotropic. 
Neglecting  the  substrate  anisotropy  can  introduce  significant  errors  in  the  design 
of  microwave  printed  circuits[3].  Similar  sensitivities  to  anisotropy  are  exhibited 
by  printed  antenna  structures.  In  some  instances  the  anisotropy  is  a  principal 
feature  of  the  design  materials,  as  is  the  case  with  biased  ferrite  substrates.  Thus, 
accurate  models,  which  include  anisotropy,  are  needed  for  microstrip  antennas  and 
arrays. 

The  analyses  for  microstrip  dipole  elements  on  general  anisotropic  substrates  have 
been  developed.  The  solutions  for  infinite  arrays  of  microstrip  dipoles  and  probe- 
fed  patches  on  generalized  anisotropic  substrates  have  also  been  derived.  Provision 
has  been  made  for  both  dielectric  and  magnetic  anisotropy.  In  Chapter  1  the 
printed  dipole  antenna  on  a  two  layer  uniaxial  substrate  is  addressed.  Chapter  2 
deals  with  dipoles  on  gyrotropic  substrates  (such  as  the  biased  ferrite).  The  use  of 
biased  ferrite  substrates  provides  very  interesting  capabilities,  including  element 
pattern  scanning,  element  pattern  shaping,  and  element  gain  enhancement.  All 
these  features  can  be  actively  controlled  by  way  of  the  bias  field  strength.  In 
addition,  in  Chapter  3  the  use  of  biased  ferrite  substrates  for  the  reduction  of  RCS 
is  illustrated.  The  non-reciprocal  nature  of  the  microstrip  antenna  with  biased 
ferrite  materials  permits  the  reduction  of  RCS  while  transmit  gain  is  preserved. 

A  substantial  effort  was  devoted  to  the  development  of  efficient  algorithms  for  the 
evaluation  of  the  spectral  integrals  that  arise  in  the  solution  method  which  was 
adopted.  It  resulted  in  significant  improvements  in  the  numerical  efficiency  of  the 
solutions.  This  is  discussed  in  Chapter  4. 

In  Chapter  5  the  analyses  for  the  infinite  phased  arrays  of  strip  dipoles  and  probe- 
fed  patches  on  general  anisotropic  substrates  are  presented.  The  results  include 
the  impedance  as  a  function  of  scan. 

The  solutions  which  were  developed  are  fullwave  integral  equation  solutions  using 
the  method  of  moments.  The  Green’s  functions  for  the  structures  were  found  by  a 
matrix  method  which  can  be  easily  extended  to  any  number  of  anisotropic  layers. 

The  effects  of  substrate  anisotropy  on  printed  antenna  characteristics  have  been 
shown  to  be  important.  In  addition,  several  potentially  significant  applications  of 
biased  ferrite  substrates  have  been  identified. 
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Chapter  1 


Microstrip  Dipole  on  a  Two-layer 
Uniaxial  Dielectric  Substrate 


The  fullwave  solution  for  the  considered  structure  involves;  (i)  the  development 
of  the  Green’s  Function  appropriate  to  the  structure;  (ii)  the  application  of  the 
method  of  moments  to  an  electric  field  integral  equation  involving  the  unknown 
currents  on  the  strip  antenna  element;  and  (iii)  the  derivation  of  the  antenna 
characteristics  such  as  input  impedance,  gain,  and  efficiency  from  the  knowledge 
of  those  currents. 


1.1  Green’s  Function 


The  geometry  of  the  structure  under  consideration  is  shown  in  Figure  1.1.  The 
bottom  layer  of  thickness  b  (region  I)  is  assumed  to  have  a  permeability  and 
permittivity  c'l 


«i  = 


( 


0 

0 


(1.1) 


0 


while  the  top  layer  of  thickness  t  (region  II)  is  assumed  to  have  a  permeability  ^2 
and  permittivity  £'2 


( 


0 


0 

0 


0 

0 


(1.2) 


The  total  material  thickness  is  h  =  (b+t).  For  the  problem  considered  here  the 
source  (antenna)  is  located  either  at  the  material  interface  (z  =  b)  or  at  the  air- 
dielectric  interfaice  (z  =  h).  Therefore,  M20cweirs  curl  equations  for  each  region 
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I 

ground  piano 


Figure  1.1:  A  Two-Layer  Printed  Antenna  Structure  with  Uniaxial  Anisotropy 


(1.3) 


are  source  free.  The  Maxwell’s  curl  equations  are 

V  X  E  =  —junH 

and 

V  X  H  =  juj€  ■  E. 

The  Fourier  and  inverse  Fourier  transform  pairs  are  defined  as: 


(1.4) 


L(x~xl,y-yl,z)  =  ~  £(A.,A„z)£>^'('-”)e'^-l''-»')<iA,<lA,  (1.5) 


and 


roo  roo 

C{Xx,Xy,z)  —  I  /  £(x  -  x/,y  -  (1.6) 

The  vector  function  £  is  either  the  electric  or  magnetic  field.  From  Eqs  1.3-  1.6, 
it  can  be  shown  that  the  transformed  field  components  miist  satisfy  the  following 
equations: 


f^-^  +  {kleg^-X^)Eg^O, 

(1.7) 

^  +  {kU,  -  X^)H.  =  0 

(1.8) 

(*oVf*  +  -~)Eg=  unXyHg  +  jA.^, 

(1.9) 

(A>.  +  ^)£,  = 

(1.10) 

(*>.  +  ^)^.  =  -“'‘.A,i.  +  I'A,^, 

(1.11) 

(*.Vf.  +  ~)H,  =  ui.X.E,  +  jX,^^‘ . 

(1.12) 

From  Eqs.  1.7-  1.12,  the  electromagnetic  fields  in  each  region  can  be  expressed 
2is: 


In  region  1,  0  <  z  <  b  , 


Eg  =  44cosh(gj2),  Hg  =  Bsinh(9jz), 

r  _  IjJHoHlXy  „  .  ^  ,  I,  ^  g  ■  \  !  a  \ 

Eg  -  — - fisinh{gjx)  +  smh(gi2), 


Ey  = 


-uHoHiX 


^5sinh(9j^)  -h  :i:^^sinh(g^) 
Eg  =  B  cosh{q[z)  -  A  cosh{qiz), 


(1.13) 

(1.14) 

(1.15) 

(1.16) 
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=  +  (1.17) 

where  - 

9?  =  \/|7(^’ -«;*>■)  ('•'«> 

,*  =  v/{A>  -  (119) 


In  region  U,  b  <  z  <  k  , 

=  Ccosh(g“z)  +  Dsinh(9“z),  H,  =  Ecosh{qlz)  +  Fsinh(g*2z),  (1.20) 

Er  =  +  Fsinh(9‘2)]  +  ^^—^[Dcoshiq^z)  +  Csinh(g“2)], 

(1.21) 

J  ^  f  sinh(^i)]  +  +  CsinhtjJz)], 

(1.22) 

=  ^^^[Fcosh(<7‘z)  +  Fsinh(9‘z)|  -  [C'cosh(g2°z)  +  Psinh(9“2)], 

^  (1.23) 

Hy  =  cosh{qlz)  +  Fsinh(922)]  +  [C  cosh(g2z)  +  Dsinh(g52)], 

^  (1.24) 


where 


In  region  0,  h  <  z  , 


where 


and 


92  =  ~  ^J*oM2) 

92  =  \/a*  -  elkln2). 


^X=( 


UHgXy  jj  ^  ^\_-Q{x-h] 


A* 


F,  =  - 


Hx  =  { 


A* 

-jA,? 

A* 


A* 

y\9 

A* 

OJCoK 


G)e- 

G-)e-9(*-'‘) 


+ 


A2 
we^A, 


G)e 


-q{t-h) 


9  =  \/^’  - 1=2 

A=yAfT^. 


(1.26) 

(1.26) 

(1.27) 

(1.28) 

(1.29) 

(1.30) 

(1.31) 

(1.32) 

(1.33) 
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For  an  x  directed  6  source  at  x  =  x/,  y  =  y/,  with  the  source  located  at  the  material 
interface,  one  has 

Hy{z  =  b^)  -  H^{z  =  b-)  = (1.34) 
If  the  source  is  located  at  the  air-dielectric  interface,  one  has 

Hy[z  =  h^)  -  Hy{z  =  h-)  = -1.  (1.35) 


Other  tangential  field  components  are  continuous  across  each  interface.  By  enforc¬ 
ing  these  boundary  conditions  and  after  some  algebraic  manipulations,  one  can 
obtain  a  system  of  linear  equations  in  the  form: 

^tanh(gifc)  — cosh{y2fc)  —  sinh{q2fc)  0 

— yjsinh(yj/>)  92Cosh(y26)  0 

0  1  tanh(95/i)  —  1/m2 

0  q2^anh(y2^)  q\  q 


B 

E 

F 

H 


=  [V'll  (1-36) 


■  -g^c5/€jtanh(qJ6) 

s\nh{q^b)qle\ltl 

cos\i{qlb)qlt\l 

0  ■ 

A  ■ 

-€^cosh(956) 

— £2sinh(q5i>) 

0 

C 

0 

tanh(y5/i)<7j£f/ef 

9 

D 

0 

tanh(g2h)£2 

-1  . 

.  G  . 

(1.37) 

For  the  case  where  the  source  is  located  at  z  =  b,  the  excitation  matrix  is 


■  0  ■ 

0 

jXy 

0 

and  = 

Ax/w£o 

0 

0 

0 

(1.38) 


While  for  the  case  where  the  source  is  located  at  z  =  h,  the  excitation  matrix  is 


■  0  ■ 

0 

0 

0 

.  -JK  . 

and[V'2j  = 

0 

0 

.  Ax/u;£<,  . 

(1.39) 


The  coefficients  of  the  Fourier  transformed  fields  can  be  found  by  inverting  the 
matrix  in  Eq.  1.36  and  1.37.  From  these  coefficients  and  Eqs.  1.5,  1.13-  1.33, 
the  electromagnetic  fields  in  the  entire  half  space  due  to  an  infinitesimal  electric 
dipole  (Green’s  function)  can  be  found. 
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1.2  The  Method  of  Moments 


Assuming  that  the  current  is  flowing  in  the  x  direction  on  the  dipole,  the  pertinent 
integral  equation  is 

=  I  ^  [  G„{x  -x\y-  y')J^{x\  y')dx'dy\  (1.40) 

J-w/2  Jo 

where  G**  is  the  Green’s  function  and  is  in  the  form; 


G.4x  -x\y-  y’)  =  ^  T  T  G„(A.,  'I  dA,.  (1.41) 

47r^  •/-oo  -00 


The  Fourier  component  G**  is  exactly  the  expression  of  Eq.  1.21  with  either 
z  =  b  (dipole  at  the  material  interface)  or  z  =  h  (dipole  at  the  air-dielectric 
interface).  In  the  method  of  moments  procedure,  the  current  on  the  dipole  is 
expanded  in  terms  of  a  set  of  overlapped  piecewise  sinusoidal  functions  and  a 
Maxwellian  function  [l].  By  assuming  that  the  dipole  length  L  is  (N  -f  l)d  and 
employing  the  method  of  moments  followed  by  a  Galerkin’s  procedure,  the  integral 
equation  is  discretized  into  a  system  of  linear  equations  which  may  be  expressed 
in  matrix  form  as  \Z]\I]  =  (V],  where  the  elements  of  [Z]  can  be  written  as 

=  I  j.  (A.)  cos(A.{m  -  n)d]d\,  ^  G„(A.,  \,)J^{\,w/2)dX.  (1.42) 


with 


cos  Ajd  —  cos  ktd 


k]  -  A» 


(1.43) 


The  impedance  matrix  elements  are  evaluated  by  a  newly  developed  method  which 
is  described  in  a  later  chapter  (6j.  This  scheme  computes  the  double  infinite 
integral  in  a  rectangular-form. 


The  elements  in  the  excitation  matrix  [V]  are 

V;  =</„,£:,  >  (1.44) 

where  n  =  1,  2,  3, .  For  a  center-fed  dipole,  supposing  that  N  is  odd,  one  has 

V+i)/2  =  1  (1-45) 

and  V*  =  0  elsewhere.  The  input  impedance  of  the  center-fed  dipole  is  defined  as 

Zin  =  V’(„+,)/2//(n+l)/2  (l-46) 
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1.3  Antenna  Gain  and  Radiation  Efficiency 


Once  the  current  in  the  dipole  and  the  Green’s  function  are  known,  it  is  possible 
to  find  the  corresponding  electromagnetic  fields  in  the  far  zone. 

The  tangential  electric  fields  at  the  air-material  interface  (z=h)  can  be  written  as 

E{x,  y,z  =  h)  =  ^rr  dK  dXy,  (1.47) 

AlC*  7—00  J -00 


where 


E{Xi,  Ay)  =  Et{Xf,  Xy)x  -I-  Ey{Xx,  Ay)j/. 


(1.48) 


j&(Ai,  Ay),  a  vector  quantity,  is  the  multiplication  of  the  spectral  dyadic  Green’s 
function  and  the  Fourier  transform  of  the  current  in  the  dipole.  The  far  zone  elec¬ 
tric  fields  are  related  to  the  tangential  spectral  electric  fields  through  the  formulae 
[5,  p.67] 


E^,  —  jko- 


Et  —  jko' 


2irR 


|.Ey  (fcl,  ky)cos<j>-  E^(k^,  fcy)sin<^]  , 
[Ez{kz,  ky)cosd>  Ey{kz,  A:y)sin(^]  , 


where 


kz  =  kosm6cos<l> 


(1.49) 

(1.50) 


(1.51) 


kz  =  kosin6sin<l> 

The  power  density  in  space  is  given  as 


=  2^  [I I’  +  I ■ 


while  the  total  radiated  power  is 


Prad  =  r'^  r  Po{0,  4>)R^  sin  0d<i>de. 
Jo  Jo 


The  transmit  gain  of  the  dipole  is  given  by 


gain(tf,  <f))  = 


P4e,  <i>)4nR^ 


and  in  Db, 

gain£,j(fl,  <f>)  =  lOlogiogain(0,  <^). 
The  radiation  efficiency  is  defined  as 


(1.52) 


(1.53) 


(1.54) 


(1.55) 


(1.56) 


(1.57) 
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where 


Ptot  =  \Rin\Iin\^  (1-58) 

For  a  lossless  material  the  total  input  power  is  the  sum  of  the  radiated  power  and 
the  power  due  to  surface  waves. 


1.4  Results 


In  this  section,  the  effects  of  material  anisotropy  on  printed  circuit  antenna  charac¬ 
teristics  are  presented  through  an  accurate  numerical  computation.  For  practl.;al 
applications,  printed  antennas  are  usually  covered  with  a  protective  layer.  This 
is  the  configuration  studied  in  this  research.  Interesting  phenomena  for  printed 
antenna  characteristics  due  to  the  presence  of  the  cover  layer  are  described  in  [2]. 
The  results  presented  here  focus  on  the  effects  of  either  substrate  or  superstrate 
anisotropy.  As  was  discussed  in  [3],  many  practical  integrated  circuit  substrates 
exhibit  strong  anisotropy  due  to  either  their  nature  or  what  is  introduced  in  the 
manufacturing  process.  Epsilam-10  (e,  =  10.2  and  c,  =  13),  Sapphire  (c,  =  9.4 
and  €t  —  11.6)  and  PTFE(CuClad)  (e,  =  2.43  and  t,  =  2.88)  are  examples  of 
these  materials. 

Typical  printed  antenna  gain  is  about  6  Db.  Figure  1.2  shows  the  broadside 
{6  =  0)  gain  as  a  function  of  frequency  with  Sapphire  as  a  substrate.  The  solid 
line  is  for  the  case  of  neglecting  the  substrate  anisotropy,  while  the  dash  line  is 
for  the  case  of  including  the  substrate  anisotropy.  It  is  seen  from  Figure  1.2 
that  the  substrate  anisotropy  has  little  effect  on  the  antenna  gain  in  almost  the 
entire  frequency  range.  On  the  other  hand,  the  superstrate  anisotropy  has  strong 
effects  on  the  gain  performance  as  illustrated  in  Figure  1.3,  where  the  anisotropic 
material  (Sapphire)  is  used  as  the  superstrate.  The  deep  valleys  occurring  in  both 
Figures  1.2  and  1.3  are  due  to  the  resonance  of  the  layered  structure.  This 
resonance  behavior  can  be  explained  in  terms  of  a  transmission  line  equivalence 
model  [2]. 

Microstrip  structures  support  surface  waves,  of  which  the  lowest  order  mode  has 
zero  cutoff  frequency.  For  printed  circuit  antenna  applications,  the  surface  wave 
power  is  considered  loss  since  all  the  power  is  intended  for  radiation  into  free  space. 
The  radiation  efficiency  characteristics  can  be  used  to  examine  the  percentage  of 
power  loss  due  to  the  surface  waves  as  compared  to  the  total  power  emanating  from 
the  antenna  element.  These  surface  waves  propagating  horizontally  are  usually 
undesirable  except  for  special  prescribed  applications.  Figure  1.4  shows  radiation 
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Figure  1.2:  Broadside  Gain  Versus  Frequency.  Antenna  Is  at  the  Material  Interface 
with  h  =  0.08  cm,  and  b  =  0.04  cm.  The  substrate  is  anisotropic  (Sapphire). 
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Figure  1.3:  Broadside  Gain  versus  Frequency.  The  antenna  is  at  the  material 
interface  with  h  =  0.08  cm,  and  b  =  0.04  cm.  The  superstrate  is  anisotropic 
(Sapphire). 
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efficiency  as  a  function  of  material  thickness  with  and  without  the  substrate  or 
superstrate  anisotropy.  It  is  seen  that  material  anisotropy  does  not  affect  the 
efficiency  much.  This  is  due  to  the  fact  that  both  surface  wave  power  and  radiation 
resistance  increase  or  decrease  together  when  the  material  constants  change. 

The  effects  of  the  substrate  anisotropy  on  the  input  impedance  of  a  center-fed 
printed  dipole  are  shown  in  Figures  1.5  and  1.6.  In  these  figures,  two  cases 
are  considered.  One  is  for  the  case  where  both  the  substrate  and  superstrate 
are  isotropic  (ci  =  2.43  and  C2  =  2.2).  The  other  is  for  the  case  of  an  isotropic 
superstrate  (62  =  2.2)  and  a  substrate  with  anisotropy  in  the  z  (vertical)  direction 
(ci*  =  2.43  and  Ci,  =  2.88).  In  addition  to  the  change  in  the  resonant  length,  the 
resonant  resistance  changes  about  10%  which  is  rather  significant.  This  implies 
that  the  neglect  of  the  material  anisotropy  can  lead  to  an  unexpected  and  severe 
input  mismatch.  It  is  also  seen  that  away  from  the  dipole  resonance,  the  material 
anisotropy  has  no  significant  effect.  The  resonant  length  of  a  center-fed  dipole 
is  also  an  important  parameter  for  the  antenna  designer.  Figure  1.7  shows  the 
dipole  resonant  length  as  a  function  of  the  layer  thickness  (or  frequency)  with  and 
without  substrate  anisotropy.  The  materials  used  in  this  illustration  are  Epsilam 
(c  =  10.2)  substrate  and  Duroid  (£  =  2.2)  superstrate.  It  is  found  that  neglecting 
the  substrate  anisotropy  results  in  an  error  in  the  resonant  length  of  more  than 
10%.  The  decrease  for  the  resonant  length  due  to  the  substrate  anisotropy  is 
because  the  material  anisotropy  in  the  z  direction  increases  the  effective  dielectric 
constant  of  the  substrate.  It  is  also  noted  that  the  increase  of  the  layer  thickness 
(or  frequency)  results  in  the  decrease  of  the  dipole  resonant  length. 

From  all  the  above  observations,  it  is  clear  that  the  material  anisotropy  has  sig¬ 
nificant  effects  on  the  near  fields  and  current  distributions  of  printed  antennas. 


1.4.1  Conclusion 

In  this  work,  we  have  developed  accurate  computer  codes  for  analyzing  two-layer 
printed  antenna  structures  capable  of  including  the  effect  of  material  anisotropy. 
The  computer  codes  include  both  near  field  and  far  field  computations.  It  hsis  been 
found  that  the  material  anisotropy  h2is  strong  effects  on  the  printed  antenna  per¬ 
formance  especially  when  the  near  fields  are  considered.  The  parameters  that  are 
affected  significantly  include  input  impedance  and  resonant  length  of  the  antennas. 
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Figure  1.4:  Radiation  Efficiency  versus  Material  Thickness.  The  antenna  is  at  the 
material  interface  with  h  =  2b. 
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Figure  1,5:  Input  Resistance  of  a  Center-fed  Dipole.  The  antenna  is  at  the  material 
interface  with  h  =  0.2  Aq,  b  =  0.1  Aq,  and  w/h  =  0.05  . 
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Figure  1.6:  Input  Reactance  of  a  Centc'-fed  Dipole.  The  antenna  is  at  the  material 
interface  with  h  =  0.2  Aq,  b  =  0.1  Aq,  and  w/h  =  0.05  . 


14 


iccie  resonant  lencth 


Figure  1.7:  Dipole  Resonant  Length  versus  Material  Thickness.  The  antenna  is  at 
the  material  interface  with  h  =  2b  and  w  =  0.05  h. 
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Chapter  2 


Microstrip  Dipoles  on  Gyrotropic 
Substrates 

The  characteristics  of  a  microstrip  dipole  on/in  a  gyrotropic  substrate  are  inves¬ 
tigated.  As  in  the  case  of  uniaxial  permittivity,  the  boundary  value  problem  is 
formulated  as  an  integral  equation  that  is  solved  by  the  method  of  moments.  The 
Green’s  function,  which  is  central  in  the  solution,  is  represented  in  the  spectral 
domain.  A  numerical  technique  for  generating  the  Green’s  function  is  described. 
The  new  algorithm  for  the  evaluation  of  the  double  infinite  spectral  integral  is 
used.  The  details  of  that  algorithm  are  presented  in  Chapter  4.  Radiation  pat¬ 
terns  (transmit  gain  patterns)  are  presented  for  the  case  of  a  dipole  on  a  biased 
ferrite  substrate.  Although  the  development  of  the  Green’s  function  is  detailed 
for  the  single  general  anisotropic  layer  substrate,  the  matrix  method  can  be  easily 
applied  to  the  problem  of  multiple  anisotropic  layers.  Computed  results  for  the 
case  of  an  isotropic  substrate  with  a  gyrotropic  (biased  ferrite)  cover  layer  are 
presented.  Those  results  illustrate  the  use  of  the  biased  ferrite  substrates  or  cover 
layers  for  element  pattern  shaping,  element  pattern  beam  steering,  and  element 
gain  enhancement. 


2.1  Introduction 

The  printed  circuit  antenna  is  an  attractive  antenna  structure  for  its  conformabil- 
ity,  low  profile,  light  weight,  and  low  cost  [7,8].  It  has  been  demonstrated  in  [1,2] 
that  the  material  constants  and  thickness  are  critical  to  the  performance  of  the 
printed  circuit  antenna.  Also,  the  frequency  range  of  operation  is  limited  (nar¬ 
row  band);  once  the  structure  is  built,  it  is  not  possible  to  change  the  antenna 
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characteristics  and  frequency  band  of  operation. 

By  using  biased  ferrite  substrates  in  microstrip  structures  one  can,  by  way  of 
the  applied  bias  dc  magnetic  field,  change  the  permeability  tensor  and  thus  the 
electrical  properties  of  the  material.  It  is  possible  to  change  and  control  the 
antenna  characteristics  dynamically  without  changing  the  physical  structure. 


2.2  Theory 

2.2.1  Material  Characteristics  of  Ferrites 


Ferrite  materials  have  been  extensively  used  in  microwave  devices  due  to  their  non¬ 
reciprocal  electrical  properties.  The  magnetic  properties  of  ferrites  arise  mainly 
from  the  magnetic  dipole  moment  associated  with  the  electron  spin.  The  funda¬ 
mental  physics  of  ferrites  have  been  well  discussed  in  the  literature,  as  for  example 
in  [9,10].  In  this  section,  we  summarize  the  material  properties  of  ferrites  in  terms 
of  a  permeability  tensor  which  is  suitable  for  the  study  of  the  wave  propagation 
problems. 


For  a  z  directed  dc  bias,  the  permeability  tensor  of  a  ferrite  material  is 


where 


Mo 


M  JK  0 

-Jk  h  0 
0  0  1 


M 


1  + 


U}QU)m 

2  2  ’ 
Wq  —  U; 


K 


UlJ  -  ’ 


Wo  = 


w™  =  , 

-y  =  1.759  X  10“''  kgfeoul , 


(2.1) 


(2.2) 

(2.3) 

(2.4) 

(2.5) 

(2.6) 


(jj  is  the  a.c.  operating  frequency,  M,  is  material  saturation  magnetization,  and  Ho 
is  the  dc  bias  magnetic  field.  For  the  geometry  shown  in  Figure  2.1,  the  planar 
ferrite  slab  is  parallel  to  the  x-y  plane.  The  dc  bieis  H  field  is,  for  practicality, 
taken  to  be  in  the  planar  direction.  For  such  cases,  the  permeability  tensor  is 


Mo 


{no-  m)  <^o 
sin  2«^o 
jKsin^o 


sin  2<f>o 

M  +  (mo  -  m)  sin*  <i>o 

-jKCOS<f>o 


-jKs'm  <i>Q 
j  K.  cos  <j)o 
M 


(2.7) 
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where  is  the  angle  of  the  direction  of  the  dc  bias  H  field  measured  from  the  x 
axis. 

It  can  be  seen  from  Eqs.  2.1  -  2.7  that  the  permeability  tensor  for  a  ferrite  is  deter¬ 
mined  by  an  inherent  material  property  (saturation  magnetization),  the  strength 
and  angle  of  the  input  bias  H  field,  as  well  as  the  operating  frequency. 


2.2.2  Green’s  Function  for  a  Grounded  Ferrite 


In  dealing  with  the  Green’s  function  for  a  grounded  ferrite,  no  more  difficulty 
is  introduced  by  considering  both  permittivity  and  permeability  as  anisotropic. 
With  this  generalization,  both  the  permittivity  and  permeability  are  tensor  forms, 
which  in  a  rectangular  Cartesian  coordinate  system  are  expressed  as: 


A  =  Mo 

Mi/z 

Mzy 

Mvv 

Mzz 

Myz 

and  c  =  to 

fzx 

CyX 

f*v 

fyy 

fzz  ' 
fyz 

Mzy 

Mzz  . 

.  G* 

Gy 

G*  . 

Then  Maxwell’s  equations  for  the  fields  inside  the  ferrite  slab  (0  <  2  <  d)  with  no 
source  are 

V  X  E  =  •  H  and  V  x  H  =  jui-  E.  (2.9) 

Introducing  the  following  Fourier  transformation 

S  =  :r  /”  a  *v) (210) 

ii  47r  J-00  li 

and  using  Eq.  2.9,  one  obtains 

-j)ixH  +  -^{zxH)=julE  (2.11) 

az 

and 

-  jk  X  E  +  — (z  X  E)  =  -jtjJfL  •  H,  (212) 

az 

with  k  =  kgX  +  kyy  and  E  and  H  are  the  spectral  electric  and  magnetic  fields, 
respectively,  more  over 

E  =  ExX  +  Eyy  +  EfZ  (213) 

and 

H  HxX  +  Hyy-\-HxZ.  (2.14) 
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After  some  algebraic  manipulations,  Eqs.  2.11  -  2.12  can  be  written  in  the  form 
of  a  set  of  linear  differential  equations,  as 


—  [^(2))  =  [A]  \^{z)] 


where 


1^(2)]  = 


and  [A]  = 


kzHx{z)  +  kyHy{z) 
kyHj^[z)  —  kgHy{z) 
kxEx{z)  +  kyEy{z) 

[  kyEx{z)  -  k^Ey\z)  J 
The  solution  of  Eq.  2.15  is 

[^(^Ol  =  [T{z\  - 
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O22 
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where 


[r(2i-22)]  =  [<^] 


gAi(ri  -*3) 

0 

0 
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0 

0 

0 


0 

0 

0 


0 

0 

0 

^>4(11 -*2) 


{<!>] 


-1 


(2.15) 


(2.16) 


(2.17) 


(2.18) 


in  which  the  A,,  i  =  1,2, 3, 4  are  the  eigenvalues  and  [^]  is  a  4  x  4  eigenvector 
matrix  of  [A].  According  to  the  Calley-Hamilton  theorem,  the  matrix  [T{zi  —  22)] 
can  also  be  expressed  as 

[7(21-22)]  =  ^  Cn{zi  -  Z2)\A]'' ,  (2.19) 

n=0 

where  the  coefficients  Cn  are  given  from  the  solution  of  the  matrix  equation 

gA,(,.-,2)  ^  Yi  \)Cu(z2  -  ^2),  j  =  1,  2,  3,  4.  (2.20) 

k=0 

The  matrix  |T]  relates  the  tangential  electromagnetic  fields  at  one  surface  (2  =  21) 
to  another  (2  =  22).  The  tangential  electromagnetic  fields  at  the  air-material 
interface  can  now  be  related  to  those  at  the  ground  plane  as 

l^.(rf-)|  =  |r(<()|  \m]  (2.21) 


The  electromagnetic  fields  in  the  air  (2  >  d)  can  be  derived  in  a  straightforward 
manner.  With  this  result  together  with  the  boundary  conditions  that  the  tangen¬ 
tial  electric  fields  are  zero  at  the  ground  plane,  one  obtains 


[0(d+)]  = 


jyjk^  -kla 

c 

cjcq  6 
r"  - 

and  1^(0)]  = 

d 

jV*!  -  ki  b 

0 

-w/io  d 

0 

(2.22) 
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where  o,  6,  c  and  d  are  unknown  spectral  quantities  to  be  determined. 


With  the  current  source  at  the  air-material  interface,  the  electromagnetic  fields 
just  above  and  below  z=d  are  related  as  follows: 

|^(rf+)|  -  =  \U,],  (2.23) 


where  [Uj]  is  either  [f/*]  or  [f/y]  and 


-ky 

kz 

If'.)  = 

kz 

0 

0 

and  \Uy\  = 

-ky 

0 

0 

(2.24) 


\Ui\  is  for  the  current  source  in  the  x  direction,  while  [t/y]  is  for  the  current  source 
in  the  y  direction.  Eq.  2.23  involves  a  4  by  4  matrix  equation,  where  the  solution 
determines  the  quantities  d,  6,  c  and  d.  As  a  result,  the  spectral  Green’s  function 
of  the  pertinent  structure  is  determined. 


2.2.3  The  Moment  Method  Solution 

The  Impedance  Matrix  Elements 

The  integral  equation  for  the  problem  is 

E,  =  j  G„  J,  ds  ,  (2.25) 

where  the  microstrip  dipole  is  oriented  in  the  x  direction,  and  is  assumed  to  have 
current  flowing  only  in  the  x  direction.  The  integral  in  Eq.  2.25  is  a  surface 
integral  over  the  strip  dipole.  G„  is  the  Green’s  function  component  which  is 
the  E,  component  due  to  a  d  current  source  in  the  x  direction.  The  spectral 
component  of  this  Green’s  function  has  been  derived  in  the  previous  section. 

In  the  moment  method  procedure,  the  current  is  expanded  in  terms  of  a  set  of 
known  basis  functions.  For  the  problem  considered  the  current  in  the  dipole  of 
length  L=  {N  +  l)a  and  width  w  is  conveniently  assumed  as 

y)  =  H  In3n{x,  y)  (2.26) 

n—\ 
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yn(x,  y)  = 


/n(x) 

jrJ{ti;/2)*  -  y* 


where 


/n(x)  = 


i;-= 


^[x  —  na|  ,  |x  —  na|  <  a 
,  otherwise. 


(2.27) 


(2.28) 


With  a  standard  moment  method  followed  by  a  Galerkin’s  procedure,  where  the 
same  expansion  and  testing  functions  are  used,  the  integral  equation  in  Eq.  2.25 
is  converted  into  a  matrix  equation  as 


In  Zmn  =  m  =  1,  2,  3  ....N, 


(2.29) 


where 


/w/2  f*i/2  rL/2  fL/2 

/  /  /  G„{x,xi,y,yi)jn{xt,yi)j„,{x,y)dxdxidydyi  (2.30) 

-iu/2  J-iu/2  J-LI2  J-L/2 


v/2  J-uj/2  J-LI2  J-LI2 


l‘iu/2  fLl2 

/  ^*i^yy)3m[x,y)dxdy.  (2.31) 

J-u)/2  J-LI2 

On  the  strip,  however.  Eg  vanishes  except  at  the  delta-gap  feed  point.  Thus,  the 
excitation  vector  terms  are  all  zero  except  that  element  that  corresponds  to  the 
testing  function  taken  over  the  delta-gap  point:. 


for  m  corresponding  to  the  feed  point 
otherwise. 


(2.32) 


The  Green’s  function  and  its  spectral  component  is  related  through  the  formula: 

G..{x,  XI,  y,  yl)  =  7^7/”/”  k,)  e-'*.!*-"-)  dk.  dk,  (2.33) 

After  inserting  Eq.  2.33  into  Eq.  2.30  and  changing  the  order  of  integration,  the 
impedance  matrix  element  Zmn  in  Eq.  2.30  may  be  written  as 

2m.  =  r  r  *,)  dk.  dk..  (2.34) 

J—oo  J—oo 


The  dipole  input  impedance  is  found  from  the  computed  current,  /,„,  at  the  delta- 
gap  feed  point: 


Zin  —  Rin  +  J  Xin  —  y- 


(2.35) 
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Far  Field  Computation 


Once  the  current  in  the  dipole  and  the  Green’s  function  are  known,  it  is  possible 
to  find  the  corresponding  electromagnetic  fields  in  the  far  zone. 

The  tangential  electric  fields  at  the  air-material  interface  (z=d)  can  be  written  as 

£(x,  dk,  dky  ,  (2.36) 

where 

E{k,,  ky)  =  ky)x  +  Eyik,,  ky)y.  (2.37) 

E{kg,  ky),  a  vector  quantity,  is  the  multiplication  of  the  spectral  dyadic  Green’s 
function  and  the  Fourier  transform  of  the  current  in  the  dipole.  The  far  zone  elec¬ 
tric  fields  are  related  to  the  tangential  spectral  electric  fields  through  the  formulae 
[5,  p  671 


E, 

=  jko  2^^  [Ey{k^,  ky)  cosd)  E^ik,,  A:Jsm<^J  , 

(2.38) 

E» 

—  jka  ky)cos<^-^-  Eyijcxi  fcy)sin(^J  , 

(2.39) 

where 

kx  =  kos\n6cos4> 

(2.40) 

and 

ky  =  kos\n6s'm<i> 

(2.41) 

The  power  density  in  space  due  to  the  current  in  the  dipole  is  given  as 

=  240.  [1  1’  ^  1  ’ 

(2.42) 

The  antenna  transmit  gain  <t>)  is 

r  rff  M  4nPtr{e,<i>) 

(2.43) 

The  receive  gain  (receiving  cross  section)  must  be  computed  from  the  results  of 
the  scattering  problem,  since  the  system  is  not  reciprocal.  The  scattering  problem 
is  addressed  in  Chapter  3  in  connection  with  the  problem  of  RCS  reduction  using 
biased  ferrite  substrates. 

The  radiation  efficiency,  r;,  for  a  lossless  substrate  is 


P,r(6 A)  dUj  ^  ^ 

iHin  |/m|’  Plot 


(2.44) 


23 


2.3  Numerical  Results 


Figures  2.2  and  2.3  are,  respectively,  El-plane  and  H-plane  transmit  gain  patterns. 
The  substrate  parameters  are  selected  to  represent  a  typical  ferrite;  the  permittiv¬ 
ity,  £r,  is  12.5,  and  the  saturation  magnetization,  moM,,  is  taken  as  0.25  Tesla.  The 
dipole  is  oriented  in  the  z  direction  which  is  parallel  to  the  <f>  =  0°  plane.  In  each 
figure  two  cases  are  plotted,  corresponding  to  two  bias  magnetic  field  directions. 
The  solid  thin  line  curves  represent  the  cases  where  the  biasing  magnetic  field  is 
y-directed,  that  is  <f>o  =  90®;  while  the  starred  line  curves  represent  the  cases  for 
which  the  field  is  z-directed,  <f)o  =  0°. 


2.3.1  Element  Pattern  Shaping,  Scanning,  and  Gain  En¬ 
hancement 


Both  for  single  element  antennas  and  antenna  arrays  the  importance  of  the  element 
factor  to  the  overall  antenna  system  performance  is  obvious.  In  many  applications 
there  would  be  a  significant  benefit  from  an  element  factor  with  higher  gain  and 
shaped  radiation  patterns.  In  scanning  arrays,  an  element  factor  with  enhanced 
off-broadside  gain  could  help  to  offset  the  loss  in  gain  that  is  associated  with  the 
decrease  in  the  projected  aperture  as  the  beam  is  scanned  off-broadside.  With 
a  sufficiently  directive  element  factor  a  single  element  could  replace  an  array  of 
elements.  If  that  single  highly  directive  element  factor  could  also  be  scanned  the 
single  element  could  replace  a  complete  phased  array. 

As  suggested  by  the  results  that  follow,  the  microstrip  antenna  element  on  a 
biased  ferrite  substrate  may  be  used  to  realize  the  single  element  characteristics 
just  described. 

Figure  2.4  is  a  plot  of  the  E-plane  directivity  pattern  of  a  strip  dipole  on  a  biased 
ferrite  substrate.  The  bias  magnetic  field  is  in  the  plane  of  the  substrate  and 
along  the  direction  of  the  dipole.  The  pattern  shows  the  enhanced  gain  at  the  off- 
broadside  angles,  with  approximately  the  same  directivity  at  60®  as  at  broadside. 
The  H-plane  pattern  is  plotted  in  Figure  2.5.  It  shows  a  striking  reshaping  of  the 
element  pattern.  The  pattern  is  highly  directive  and  the  element  pattern  “beam” 
goes  out  at  a  grazing  angle.  Those  patterns  were  computed  at  30  GHz  and  for  a 
fixed  bias  magnetic  field  of  7150  Gauss.  The  substrate  was  2.0  mm  thick  and  with 
a  relative  permittivity  of  12.6.  The  saturation  magnetization  of  the  material  was 
3000  Gauss. 
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Figure  2.2:  E-plane  Transmit  Gain  Pattern  for  a  Center-fed  Dipole  on  a  Ferrite 
Substrate. 

The  dipole  is  oriented  in  the  <f>o  =  0®  direction.  Substrate  thickness  =  0.04  A^, 
dipole  width  =  0.02  A,,,  =  12.5,  n  =  0.766,  and  k  =  -0.567  . 

Two  cases  are  plotted: 

(1)  ****  the  bias  field  angle  is  </»<,  =  0®,  the  dipole  resonant  length  is  =  .366Ao  ; 

(2)  —  the  bias  angle  is  90°,  the  dipole  resonant  length  is  1,  =  .256Ao  . 
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Figure  2.3:  H-plane  Transmit  Gain  Pattern  for  a  Center-fed  Dipole  on  a  Ferrite 
Substrate. 

The  dipole  is  oriented  in  the  (f>o  =  0°  direction.  Substrate  thickness  =  0.04  Xg, 
dipole  width  =  0.02  A,,,  =  12.5,  //  =  0.766,  and  k  =  -0.567  . 

Two  cases  are  plotted: 

(1)  ****  the  bias  field  angle  is  <po  =  0°,  the  dipole  resonant  length  is  U  =  .366Ao  ; 

(2)  —  the  bias  angle  is  90°,  the  dipole  resonant  length  is  /,.  =  .256Ao  . 


26 


20  10  0  -10  -20  -30  -20  -10  0  10  20  dB 


Figure  2.4:  E-plane  Directivity  for  a  Center-fed  Dipole  on  a  Ferrite  Substrate. 

The  frequency  is  30  GHz.  The  bias  field  is  parallel  to  the  dipole. 

The  substrate  thickness  is  2.0  mm,  and  =  12.6. 

Hde  =  7150  Gauss,  and  M,  =  3000  Gauss. 
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Figure  2.5:  H-plane  Directivity  for  a  Center-fed  Dipole  on  a  Ferrite  Substrate. 

The  frequency  is  30  GHz.  The  bias  field  is  parallel  to  the  dipole. 

The  substrate  thickness  is  2.0  mm,  and  fr  =  12.6. 

Hde  =  7150  Gauss,  and  M,  =  3000  Gauss. 
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By  varying  the  applied  bias  magnetic  field  strength  one  can  vary  the  patterns  of  a 
printed  antenna  on  a  ferrite  substrate.  For  the  case  where  the  bias  field  is  in  the 
direction  of  the  dipole,  the  H-plane  element  pattern  “beam"  will  scan  as  a  function 
of  the  applied  bias  field.  This  is  illustrated  in  the  plot  of  Figure  2.6,  for  which 
the  configuration  consists  of  a  two  layer  structure  with  an  isotropic  substrate  and 
a  biased  ferrite  cover  layer.  The  isotropic  substrate  is  1  mm  thick  with  a  relative 
permittivity  of  10.0,  and  a  permeability  of  1.0.  The  ferrite  cover  layer  is  3  mm 
thick  with  a  permittivity  of  12.6  and  a  saturation  magnetization  (  )  of  2-500 

Gauss.  The  bias  magnetic  field  is  directed  parallel  to  the  strip  dipole.  Figure  2.6 
is  for  the  fixed  frequency  of  9  GHz.  It  is  a  plot  of  the  element  beam  pointing  angle 
as  a  function  of  the  bias  field  strength.  The  H-plane  patterns  for  two  different  bias 
field  strengths  are  shown  in  Figure  2.7.  Each  pattern  corresponds  to  one  point 
in  the  curve  of  Figure  2.6.  Figure  2.8  is  a  plot  of  the  beam  pointing  angle  as  a 
function  of  frequency  for  a  fixed  bias  field  strength  H^c  —  1000  Gauss.  From  the 
preliminary  computations  which  have  been  made  it  appears  that  the  nonreciprocal 
properties  of  the  biased  ferrite  are  more  apparent  when  the  ferrite  is  used  as  a  cover 
layer. 


2.3.2  Conclusion 

Analyses  for  microstrip  dipole  elements  on/in  gyrotropic  (biased  ferrite)  substrates 
have  been  developed  using  a  rigorous  moment  method  solution.  The  applied  bias 
field  can  have  a  significant  effect  on  the  electrical  performance.  The  applied  bias 
can  be  used  to  dynamically  change  the  dipole  characteristics,  including  resonant 
frequency  and  pattern  shape.  A  highly  directive  element  pattern  can  be  achieved. 
Moreover,  the  bias  magnetic  field  can  be  used  as  a  control  of  the  element  fac¬ 
tor  beam  direction  (peak).  This  suggests  that  the  structure  can  be  used  as  an 
electronically  tunable  (adaptable)  multifunctional  system. 
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Figure  2.6:  H-plane  element  beam  pointing  angle  as  a  function  of  the  bias  magnetic 
field  strength. 


The  frequency  is  9  GHz.  The  bias  field  is  parallel  to  the  dipole. 
Substrate:  thickness  =1.0  mm,  e,  =  10.0,  and  Hr  =  1.0. 


■^dc  —  1000  Oa.uss 
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Figure  2.7:  H-plane  patterns  for  H^c  =  600G  and  H4e=  1000  G. 

The  frequency  is  9  GHz.  The  bias  field  is  parallel  to  the  dipole. 
Substrate:  thickness  =1.0  mm,  =  10.0,  and  Hr  =  1-0. 

Ferrite  Cover  Layer:  thickness  =  3.0  mm,  tr  =  12.6,  M,=  2,500  G. 
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Figure  2.8:  H-plane  element  beam  pointing  angle  as  a  function  of  frequency. 


The  bias  field  strength  is  Hde=  1000  G.  The  bias  field  is  parallel  to  the  dipole. 
Substrate:  thickness  =1.0  mm,  fr  =  10.0,  and  Hr  =  1.0. 

Ferrite  Cover  Layer:  thickness  =  3.0  mm,  =  12.6,  Mt=  2500  G. 


Chapter  3 


Scattering  and  RCS  of  Microstrip 
Dipoles  on/in  Biased  Ferrite 
Substrates 


In  many  applications  the  antenna  scattering  or  radar  cross  section  (RCS)  is  of 
importance.  In  low  observable  applications  it  can  become  the  most  important 
performance  parameter.  Designs  with  simultaneous  low  in-band  RCS  and  high 
gain  and  efficiency  may  be  required.  Conventional  antenna  designs  are  limited  in 
the  extent  to  which  the  EM  scattering  can  be  reduced  by  their  reciprocal  charac¬ 
teristics.  Since  because  of  the  system  reciprocity,  a  reduction  in  antenna  in-band 
RCS  implies  a  reduction  in  antenna  gain.  However,  the  use  of  non-reciprocal  ma¬ 
terials,  such  as  biased  ferrites,  may  provide  substantial  RCS  reduction  while  the 
antenna  gain  is  preserved. 

The  antenna  scattering  problem  is  of  interest  for  a  second  reason:  it  provides  the 
means  by  which  the  receive  function  of  the  antenna  (the  receive  cross  section) 
can  be  characterized.  In  non- reciprocal  systems  the  transmit  and  receive  patterns 
are  distinct,  and  must  be  analyzed  separately.  The  characterization  of  the  receive 
antenna  involves  the  solution  of  the  problem  of  scattering  from  a  loaded  antenna. 
Specifically,  it  requires  one  to  find  the  induced  currents  on  the  loaded  antenna  when 
a  plane  wave  is  incident  on  the  structure,  which  is  the  essence  of  the  scattering 
problem. 

In  this  chapter  the  analyses  and  methods  described  in  the  previous  chapters  are 
applied  to  the  scattering  problem.  The  computation  of  RCS  for  the  microstrip 
antenna  on  a  layered  structure  including  gyrotropic  layers  is  outlined. 
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3.1  The  Integral  Equation  for  Scattering 


The  integral  equation  for  the  problem  is 

+  j  Gx,  Jx  ds  =  0,  at  the  strip  dipole,  (3.1) 

where  the  microstrip  dipole  is  oriented  in  the  x  direction  and  is  assumed  only  with 
current  flowing  in  the  x  direction.  The  integral  in  Eq.  3.1  is  a  surface  integral 
over  the  strip  dipole.  Gxx  is  the  Green’s  function  component  which  is  the  E* 
component  due  to  a  6  current  source  in  the  x  direction.  The  spectral  component 
of  this  Green’s  function  has  been  derived  in  the  previous  chapters.  E’”'  is  the  x 
component  of  the  total  electric  field  which  results  when  a  plane  wave  is  incident  on 
the  layered  structure  without  the  strip  dipole.  This  term  thus  includes  reflected 
or  scattered  plane  waves. 


3.2  The  Moment  Method  Solution  for  the  In¬ 
duced  Dipole  Currents 


In  the  moment  method  procedure,  the  current  is  expanded  in  terms  of  a  set  of 
known  basis  functions.  For  the  problem  considered  the  current  on  the  dipole  of 
length  L  =  {N  +  l)a  and  width  w  is  assumed  to  have  the  form 


N 


n=l 


where 


and 


y)  =  Z]  ^nJnix,  y) 

fn{x) 


Jn(x,  y]  = 


7r\/(u;/2)2  -  y*’ 


/n(x)  = 


(:-• 


—  ijx  —  no|  ,  |x  —  na|  <  a 
,  otherwise. 


(3,2) 


(3.3) 


(3.4) 


In  Galerkin’s  procedure  the  testing  functions  are  chosen  to  be  the  same  as  the 
expansion  functions.  The  integral  equation  (Eq.  3.1)  is  converted  into  a  matrix 
equation: 


N 

I 

n=l 


£  In  Zmn  =  ,  m  =  1,  2,  3  ....iV, 


(3.5) 


where 


fniji  /•tu/2  /•L/2  rLji 

Zmn=  /  /  /  G„{x,xl,y,yl)j„{xt,yt)jm{x,y)dxdxldydyl  {3.6) 

J~w/2  J-v)/2  J-LI2  J-L/2 
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and 


V„^  =  f  f  y)j^{x,  y)  dxdy.  (3.7) 

J-w/2  J-LI2 

The  Green’s  function  and  its  spectral  component  is  related  through  the  formula: 

G«(x,  X/,  y,  yi)  =  ^  k,)  «-»*«('-")  e->*.(v-vO  dk,  dky  (3.8) 

After  inserting  Eq.  3.8  into  Eq.  3.6  and  changing  the  order  of  integration,  the 
impedance  matrix  element  Z^n  in  Eq.  3.6  may  be  written  as 

Zmn  =  r  r  G\4k,,  J^{kywl2)  (3.9) 

^-00  J—oo  K* 


3.3  Scattered  Far  Fields  and  Radar  Cross  Sec¬ 
tion 


Once  the  current  on  the  dipole  and  the  Green’s  function  are  known,  it  is  possible 
to  find  the  corresponding  electromagnetic  fields  in  the  far  zone. 

The  tangential  electric  fields  at  the  air-material  interface  (z=d)  can  be  written  as 

E{x,  y,z  =  d)  =  ^r  r  E{k,,  dk,  dky,  (3.10) 

47r*  J -oo  J -00 

where 

E{kz,  ky)  -  E^{kr,  ky)x  +  Ey(fc*,  ky)y.  (3.11) 

E{kx,  ky),  a  vector  quantity,  is  the  multiplication  of  the  spectral  dyadic  Green’s 
function  and  the  Fourier  transform  of  the  current  in  the  dipole.  The  far  zone  elec¬ 
tric  fields  are  related  to  the  tangential  spectral  electric  fields  through  the  formulae 

15] 

g-jkoR  _ 

E^  =  jko  ky)  cos  <l>  —  Eg{^kg,  ^y)sin<^j  ,  (3.12) 

.  .  . 

*'  )cos<^-|- Ey(fcj,  /i:y)sin^J  ,  (3.13) 

where 


kg  =  A:osin9cos<^ 


(3.14) 


and 


ky  =  ilcosin^sin^ 


The  power  density  in  space  due  to  the  current  in  the  dipole  is  given  as 


(3.15) 


(3.16) 
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while  the  power  density  of  the  incident  plane  wave  is  l/(2407r). 
radar  cross  section  (RCS)  in  square  meters  is 

RCS(«,  <k)  =  4ir  [|  £,  I’  +  I  Et  I’]  , 

and  in  dBsm, 

RCSdB  =  10log,oRCS(d,  <i>). 


Therefore,  the 

(3.17) 

(3.18) 


3.4  Results 


As  was  noted  previously,  the  objective  in  RCS  reduction  schemes  is  to  reduce  the 
RCS  but  at  the  same  time  preserve  the  antenna  RF  performance,  in  particular  to 
preserve  transmit  gain  and/or  receive  cross  section. 

From  the  data  that  follows  one  can  conclude  that  with  biased  ferrite  cover  layers 
it  may  be  possible  to  reduce  RCS  and  at  the  very  least  preserve  transmit  gain. 
(The  interplay  of  RCS  reduction,  transmit  gain,  and  receive  cross  section  has  yet 
to  be  fully  investigated.) 

The  developed  solution  is  based  on  a  fullwave  method  that  takes  complete  account 
of  the  surface  wave  phenomena  due  to  the  substrate  construction.  Consequently, 
in  the  results  that  follow  the  computed  antenna  gains  include  the  surface  wave 
efficiency  factors. 

Figure  3.1  consists  of  plots  of  the  computed  RCS,  at  broadside,  for  a  microstrip 
dipole  on  an  isotropic  substrate  with  and  without  a  biased  ferrite  cover  layer.  For 
each  of  the  two  cases  the  dipole  was  terminated  with  an  appropriate  matched  load 
at  15  GHz.  The  isotropic  substrate  was  0.72  mm  thick,  with  a  permittivity  of  2.2. 
The  ferrite  cover  layer  was  1.08  mm  thick,  with  an  assumed  permittivity  of  12.6, 
and  a  saturation  magnetization  (  M,  )  of  2780  Gauss.  The  bias  field  was  2500 
Gauss  and  it  was  in  a  direction  normal  to  the  layer  surface.  The  dashed  curve, 
which  corresponds  to  the  case  with  the  biased  ferrite  cover,  shows  a  reduction  of 
the  peak  RCS  of  about  7  dBsm  relative  to  the  uncovered  isotropic  case. 

The  reduction  in  RCS  is  significant  only  if  it  is  possible  to  simultaneously  maintain 
the  gain  of  the  antenna.  This  is  in  fact  the  case.  Figure  3.2  consists  of  plots  of  the 
computed  transmit  gains  for  the  structures  with  and  without  the  biased  ferrite 
cover.  It  is  seen  that  at  15  GHz  the  gain  of  the  antenna  with  the  biased  ferrite 
cover  only  dropped  by  1.0  dB  compared  to  the  uncovered  isotropic  case.  The 
biased  ferrite  cover  layer  reduced  the  peak  antenna  RCS  by  7  dBsm  while  only 
reducing  the  transmit  gain  by  less  than  1.0  dB.  In  some  applications  that  would 
be  an  excellent  trade-off. 
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In  Figure  3.3  the  computed  broadside  RCS  is  plotted  for  the  full  range  of  incident 
field  polarizations.  This  is  for  the  same  parameters  as  those  of  Figure  3.1,  at  15 
GHz.  From  this  it  is  clear  that  the  reduction  in  RCS  is  not  just  due  to  a  rotation 
of  the  incident  field  relative  to  the  polarization  of  the  dipole.  The  worst-case 
(largest)  RCS  for  the  antenna  with  a  biased  ferrite  cover  is  lower  by  6dBsm  than 
the  worst-case  RCS  of  the  isotropic  antenna. 

Figures  3.4  and  3.5  are,  respectively,  plots  of  the  current  at  the  load  and  the 
received  power  as  a  function  of  the  incident  field  polarization.  Those  are,  as  in 
the  previous  cases,  for  broaside  incidence  at  15  GHz. 

The  curves  of  Figures  3.6  and  3.7  are  of  antenna  RCS  over  frequency  and  for 
several  different  incident  field  polarizations.  The  corresponding  induced  dipole 
currents  are  plotted  in  Figures  3.8  and  3.9. 

Computations  were  made  for  several  discrete  angles  of  incidence  other  than  broad¬ 
side.  The  data  of  Table  3.1  indicates  that  the  RCS  reduction  noted  in  the  previous 
data  is  not  specific  to  the  broadside  incidence  condition,  but  rather  that  the  phe¬ 
nomenon  is  observed  at  all  angles  of  incidence. 

Transmit  gain  and  receive  gain  patterns  are  provided  in  Figures  3.10  and  3.11, 
respectively.  The  sum  of  the  data  suggests  that  the  reduction  in  RCS  is  obtained 
at  the  expense  of  the  receive  gain  (receive-cross-section)  of  the  antenna.  This  at 
least  is  the  case  for  the  construction  and  parameters  which  were  used  here.  It  may 
be  possible  to  find  a  set  of  parameters  for  which  the  receive  gain  is  preserved  and 
instead  the  transmit  gain  is  sacrificed  in  order  to  achieve  the  reduction  of  antenna 
RCS. 


3.5  Conclusion 


The  scattering  and  RCS  for  microstrip  antennas  on  gyrotropic  substrates  has  been 
investigated.  A  comparison  of  gain  and  RCS  data  for  microstrip  antennas  with 
biased  ferrite  cover  layers  indicates  that  it  is  possible  to  reduce  broadside  RCS 
while  preserving  transmit  gain.  These  are  preliminary  results  and  a  variety  of 
parametric  studies  remain  to  be  carried  out.  The  interplay  of  RCS,  transmit  gain, 
and  receive  cross  section  must  be  fully  explored. 
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No  cover  layer,  L  =  0.706  cm  and  w  =  0.04  cm 

With  a  ferrite  cover  layer,  L  =  0.343  cm  and  w  =  0.04  cm 


10  11  12  13  14  15  16  17  18  19  20 

Frequency  in  GHz 


Figure  3.1:  RCS  at  Broadside  as  a  Function  of  Frequency  With  and  Without  a 
Biased  Ferrite  Cover  Layer. 

Substrate:  thickness  =  0.72  mm,  and  t,  =  2.2. 

Ferrite  Cover  Layer:  thicknes3=  1.08  mm,  Cr  =  12.6, 

Saturation  magnetization  (  M,  )  =  2780  G, 

=  2500  Gauss,  and  the  bias  field  is  normal  to  the  layers  and  dipole. 

In  each  case  the  dipole  is  terminated  with  a  matched  load  at  15  GHz. 


Figure  3.2:  Transmit  Gain  at  Broadside  as  a  Function  of  Frequency  for  a  Microstrip 
Dipole  With  and  Without  a  Biased  Ferrite  Cover  Layer. 

Substrate:  thickness  =  0.72  mm,  and  Cr  =  2.2. 

Ferrite  Cover  Layer:  thickness=  1.08  mm,  Cr  =  12.6, 

Saturation  magnetization  (  M ,  )  =  2780  G, 

Hdc  =  2500  Gauss,  and  the  bias  field  is  normal  to  the  layers. 


E  field  polarized  angle  in  degrees 


Figure  3.3:  RCS  at  Broadside  as  a  Function  of  Incident  Field  Polarization:  With 
and  Without  a  Biased  Ferrite  Cover  Layer. 

Frequency=  15  GHz 

Substrate:  thickness  =  0.72  mm,  and  (r  =  2.2. 

Ferrite  Cover  Layer:  thickness=  1.08  mm,  e,  =  12.6, 

Saturation  magnetization  (  M,  )  =  2780  G, 

=  2500  Gauss, 

The  bias  field  is  vertical  (normal  to  substrate). 

The  dipoles  are  terminated  with  matched  loads: 

Z/"  =  13.0n,  and  =  2.800. 
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Figure  3,4:  Current  at  the  Load  as  a  Function  of  the  Incident  Field  Polarization, 
for  a  Broadside  Incident  Field:  With  and  Without  a  Biased  Fer’^ite  Cover  Layer. 

Frequency=  15  GHz 

Substrate:  thickness  =  0.72  mm,  and  =  2.2. 

Ferrite  Cover  Layer:  thickness=  1.08  mm,  c,  =  12.6, 

Saturation  magnetization  [  M,)  —  2780  G, 

Hde  =  2500  Gauss, 

The  bias  field  is  vertical  (normal  to  substrate), 

The  dipoles  are  terminated  with  matched  loads: 

Z/"  =  13.0n,  and  =  2.800. 
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Figure  3.5:  Received  Power  as  a  Function  of  the  Incident  Field  Polarization,  for  a 
Broadside  Incident  Field:  With  and  Without  a  Biased  Ferrite  Cover  Layer. 

Frequency=  15  GHz 

Substrate;  thickness  =  0.72  mm,  and  €r  =  2.2. 

Ferrite  Cover  Layer:  thickness=  1.08  mm,  c,  =  12.6, 

Saturation  magnetization  {  M,  )  =  2780  G, 

Hde  =  2500  Gauss, 

The  bias  field  is  vertical  (normal  to  substrate). 

The  dipoles  are  terminated  with  matched  loads: 

Z/*’’  =  13.0n,  and  =  2.80n. 
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RCS  in  dBsm 


Figure  3.6:  Broadside  RCS  as  a  Function  of  Frequency  for  Different  Incident  Field 
Polarizations.  The  angle  4>  identifies  the  electric  field  direction  in  the  x-y  plane. 
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Fieciueiicy  in  GFIz 


Figure  3.9:  Induced  Current  as  a  Function  of  Frequency  for  Different  Incident 
Field  Polarizations.  The  angle  4>  identifies  the  electric  field  direction  in  the  x-y 
plane. 
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Incident  Angle 

Dielectric 

Ferrite/Dielectric 

e  =  30°,  (f>  =  90° 

-33.6 

-39.3 

e  =  60°,  <t>  =  90° 

-43.1 

-44.8 

e  =  30°,  (/i  =  0° 

<  -90 

-43.8 

e  =  60°,  =  0° 

<  -90 

-58.9 

TE  MODE  INCIDENT 


Incident  Angle 

Dielectric 

Ferrite /Dielectric 

e  =  30°,  0  =  90° 

<  -90 

-40.4 

d  =  60°,  <f)  =  90° 

<  -90 

-43.1 

d  =  30°,  0=0° 

-33.9 

-40.1 

e  =  60°,  0=0° 

-40.6 

-54.6 

TM  MODE  INCIDENT 


Table  3.1:  RCS  for  Off-broaside  Angles  of  Incidence  and  TE  and  TM  Incident 
Field  Polarizations. 
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f'  =  14  5  GHz  ^  -  15  GHz 


Figure  3.10:  Transmit  Gain  Patterns:  0  =  0"  plane. 
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Chapter  4 


Evaluation  of  Spectral  Integrals 
in  the  Moment  Method 


A  numerical  algorithm  for  the  moment  method  solution  of  printed  circuit  antenna 
problems  is  developed.  The  formulation  employs  a  spectral  domain  approach; 
however,  the  integration  variables  are  in  the  cartesian  instead  of  the  polar  co¬ 
ordinate  system.  This  numerical  scheme  separates  the  integration  involving  the 
longitudinal  and  the  transverse  dependence  of  the  basis  functions.  This  aspect 
may  reduce  the  computation  effort  drastically.  The  new  approach  is  found  very 
efficient  and  easy  to  use,  and  is  suitable  for  a  variety  of  basis  functions  including 
mixed  type  basis  function  sets. 


4.1  Introduction 


Method  of  moments  solutions  of  integral  equations  have  been  applied  extensively 
to  printed  circuit  antenna  problems  in  the  past.  In  this  approach,  an  exact  Green’s 
function  for  a  grounded  dielectric  substrate  is  used.  The  method  of  moments  trans¬ 
forms  the  integral  equation  into  linear  simultaneous  equations  where  the  solution 
(current  on  the  microstrip)  can  be  found  by  matrix  inversion.  The  efficiency  and 
accuracy  of  the  solution  depends  on  the  computation  of  each  individual  matrix 
element  as  well  as  the  matrix  inversion.  Depending  on  the  method  of  formulation, 
the  computational  methods  thus  lead  to  two  different  branches.  If  the  formulation 
is  in  a  real  space,  the  computation  involves  the  evaluation  of  a  Sommerfeld-type 
integral  [11,1,12,13].  Numerical  techniques  for  efficient  computation  of  this  type 
of  integral  are  developed  in  [14,15,16].  These  numerical  methods  have  a  common 
feature  in  that,  if  a  true  Galerkin  procedure  is  used,  then  the  computations  of  a 
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Sommerfeld  type  integral  (Green’s  function)  and  four  finite  integrals  (two  surface 
integrals  over  the  planar  conductor)  are  required.  An  efficient  approach  as  was 
described  in  [14]  is  to  tabulate  the  Green’s  function  as  a  function  of  the  distance 
between  the  source  and  field  point  such  that  the  infinite  integral  needs  to  be  com¬ 
puted  only  once.  The  four  finite  integrals  can  usually  be  reduced  two  integral  if 
the  expansion  and  testing  functions  are  properly  chosen.  For  irregular  shaped  or 
entire  domain  basis  functions,  more  analytic  derivation  may  be  required  to  reduced 
the  computation  effort  of  the  four  integrals. 

If  the  formulation  of  the  problem  is  in  a  plane  wave  spectral  domain,  the  compu¬ 
tation  always  involves  a  double  infinite  integration  which  is  usually  transformed 
into  a  finite  and  an  infinite  one  by  a  transformation  from  a  cartesian  to  a  polar 
coordinate  system  as  described  in  [17|-[23|.  The  problem  with  this  computation  is 
that  the  finite  integration  involves  an  oscillatory  function,  while  the  infinite  one  is 
slowly  convergent  and  may  be  highly  oscillatory  when  the  basis  function  spacing  is 
large.  These  aspects  lead  to  a  numerical  computation  which  is  very  inefficient.  An 
improved  technique  using  an  equivalent  image  extraction  method  has  been  pro¬ 
posed  [24].  Although  it  enhances  the  computational  efficiency,  it  is  only  developed 
for  piecewise  sinusoidal  basis  functions  with  no  transverse  variation  and  still  has 
difficulties  for  the  thin  substrate  cases. 

In  this  communication,  a  new  algorithm  for  the  matrix  element  computation  in  the 
method  of  moments  formulation  is  presented.  In  the  new  approach,  the  integration 
is  carried  out  directly  in  a  plane  wave  spectral  domain,  which  involves  a  double 
infinite  integration.  The  advantages  of  computing  the  double  infinite  integration 
directly  rather  than  treating  a  finite  and  an  infinite  one  will  be  discussed.  The 
comparison  of  the  new  and  the  previous  methods,  in  terms  of  computational  effi¬ 
ciency,  is  also  made  through  the  examples  of  a  center-fed  dipole  and  a  microstrip 
open-end  discontinuity  on  an  isotropic  substrate. 


4.2  The  New  Algorithm 


If  developed  in  a  spectral  domain  [l7]-[24],  the  impedance  matrix  elements,  from 
the  method  of  moments  formulation  of  the  printed  circuit  structure  shown  in 
Figure  4.1,  are  in  the  form  of 

z  =  r  r  (4.1) 

J  —00  j  -‘OO 

where  (xm,yn)  is  the  center-to  center  displacement  between  the  expansion  and 
testing  functions.  The  Green’s  function  G(Ajj,Ay)  includes  the  information  of  the 
geometry  and  material  constants.  If  a  single  layer  microstrip  structure  is  used  as 
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Figure  4.1;  A  Two-layer  Printed  Antenna  Structure 
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an  example,  the  Green’s  function  may  be  any  of  the  following  three  components, 
depending  on  which  current  and  electric  field  components  are  considered  [25]. 


/m  /\  \  \  ■^*1  ^r)  1 

~  '  C,(A)  +  A(A)Z)„(A)I’ 

(4.2) 

_  AjAy  A„Axgi(l  —  £r)  1 

‘A(A)  2>e(A)P«(A) 

(4.3) 

,€r-Aj  A59l(l-Cr), 

'd,(A)  D,(A)I>„(A)*’ 

(4.4) 

where 

A  (A)  -  9  +  tanh^27rgih) 

(4.5) 

D™(A)  - 

(4.6) 

II 

1 

I-* 

(4.7) 

9l  =  \/A*  -  tr 

(4.8) 

and 

A  =  v/AJ  +  AJ. 

(4.9) 

In  this  work,  we 

use  an  approach,  in  which  the  integral  in  Eq. 

4.1  is  evaluated 

directly  instead  of  transforming  it  into  a  polar  coordinate  system, 
ment,  Eq.  4.1  can  be  written  as 

With  rearrange- 

Z  =  r  Ja(Ax)c'^*”*5(Ax)dAx 

J-oo 

(4.10) 

where 

fPO 

5(Ax)  =  /  G(Ax,A„)Jt(A„)c^''»«'’*dA,. 

(4.11) 

Both  the  Ax  and  A„  integration  ranges  can  usually  be  reduced  to  the  semi-infinite 
ranges  from  0  to  oo  depending  on  the  even  or  odd  symmetry  of  the  integrand. 
The  Green’s  functions  contain  singularities  in  the  range  1  <  A  <  yJTr  which 
correspond  to  either  TM  surface  wave  modes  (i?„,(A)  =  0)  or  TE  surface  wave 
modes  {I?e(A)  =  0)  [llj.  A  pole  extraction  method  was  used  [11,1,19,20,21]  in 
the  past,  which  in  addition  to  the  numerical  integration,  requires  the  calculation 
of  residues  and  Cauchy  principal  values  at  the  singularity  points.  For  a  multi¬ 
layer  or  anisotropic  structure,  the  location  of  the  singularities  and  the  derivative 
of  Dt{X)  and  D^i^)  usually  introduces  more  complexity  if  not  more  difficulty 
into  the  problems.  A  nice  way  to  avoid  this  singularity  problem  is  to  deform  the 
integration  in  the  range  from  0  to  [22].  For  the  present  approach,  a  careful 
study  shows  that  both  the  A*  and  A^  integration  contour  needs  to  be  deformed 
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Im(A,  or  A„) 


Figure  4.2:  Integration  Contour  in  either  the  Aj  or  A^  complex  plane. 


or  \) 
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to  avoid  the  singularities.  The  deformed  integration  path  is  shown  in  Figure  4.2. 
Although  the  deformed  integration  path  can  be  arbitrary  mathematically,  care 
must  be  taken  in  choosing  the  proper  path  for  numerical  integration.  When  the 
path  is  off  the  real  axis,  the  integrand  increases  exponentially;  therefore  if  the 
contour  is  too  far  from  the  real  axis,  numerical  instability  occurs.  Also  if  the 
path  is  too  close  to  the  real  axis  (the  singularities),  the  integrand  is  not  a  smooth 
function.  A  numerical  study  for  the  integration  path  in  Figure  4.2  shows  that 
choosing  the  maximum  of  Im(A*)  and  Im(Ay)  to  be  about  l/(50h)  (h  is  substrate 
thickness)  is  adequate. 

For  A^  >  the  function  5(A,)  is  an  integral  where  the  integrand  is  a  smooth 
function  unless  y„  is  large.  Also,  the  integral  is  uniformly  convergent.  For  a  sin¬ 
gle  layer  microstrip  structure  with  a  roof-top  basis  function,  when  A  :»  the 
Green’s  function  approaches  its  asymptotic  form  which  contains  only  algebraic 
functions  and  then  integral  can  be  evaluated  either  in  a  closed  form  or  can  be 
deformed  to  an  integral  with  the  integrand  decaying  exponentially.  An  example 
of  this  derivation  is  given  in  Appendix  4-A.  If  the  basis  functions  are  not  simple 
trigonometric  functions  (for  example,  Maxwellian  function,  irregular  shape  ba¬ 
sis  functions)  or  the  structure  is  multi-layered  or  anisotropic,  there  is  no  simple 
asymptotic  form  for  the  integrand  in  5(Aa).  In  these  cases,  a  Filon’s  integration 
scheme  (26)  may  be  applied,  where  the  integration  path  is  divided  into  many  inter¬ 
val  and  the  non-sinusoidal  term  in  the  integrand  is  approximated  parabolically  in 
each  interval.  Once  the  parabolic  function  approximation  is  used,  the  integration 
in  each  divided  interval  can  be  evaluated  in  a  closed  form. 

One  distinct  feature  of  the  new  algorithm  is  that  the  function  ■S(Ai)  is  a  smooth 
and  monotonic  function  when  A*  >  y/e^.  From  Eqs.  4.2-  4.9,  it  is  seen  that 
the  Green’s  function  is  asymptotically  linear  with  respect  to  A*.  This  implies 
that  one  can  use  only  a  few  sampling  points  to  tabulated  the  function  5(A,) 
and  that  interpolation  can  provide  the  value  with  excellent  accuracy.  A  typical 
example  of  the  functional  behavior  of  5(A,)  is  shown  in  Figure  4.3,  where  it 
is  seen  that,  asymptotically,  5(A*)  is  almost  a  straight  line.  Once  the  function 
5(Ax)  is  tabulated,  one  can  proceed  to  compute  the  impedance  matrix  elements 
shown  in  Eq.  4,10,  which  contain  a  smooth  function  multiplying  either  cos(AxXto) 
or  sin(AxXTO).  This  sinusoidal  function  may  oscillate  rapidly  when  Xm  is  large. 
To  alleviate  this  difficulty,  one  can  again  approximate  the  smooth  function  with 
piecewise  parabolic  functions. 

The  appro2u:h  proposed  here  has  a  feature  that  for  fixed  the  term  <S'(Ai)  needs 
to  be  computed  and  tabulated  only  once  and  is  stored  in  an  array.  For  example, 
for  a  rectangular  patch  with  one  dimensional  current,  which  is  divided  into  M  x  N 
subsections,  using  the  new  algorithm  to  gene!  ate  the  impedance  matrix  elements, 
it  is  required  to  compute  a  one  dimensional  infinite  integration  only  M-I-NP  times 
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where  P  is  the  number  of  sampling  points  for  5{Ai).  In  contrast,  if  the  integration 
is  performed  in  the  polar  coordinate,  it  is  required  to  compute  a  double  integration 
M  X  N  times.  This  feature  enhances  drastically  the  computational  eflSciency.  It  is 
also  noted  that  the  approach  proposed  in  this  communication  may  not  require  the 
use  of  the  asymptotic  form  of  the  Green’s  function  unless  for  some  special  cases. 
In  other  words,  the  efficiency  of  the  numerical  integration  does  not  depend  on  how 
far  the  integration  path  goes  for  the  integrand  to  approach  its  asymptotic  form. 
Therefore,  the  algorithm  is  efficient  over  almost  the  entire  frequency  range. 


4.3  Numerical  Examples 

4.3.1  A  Center-fed  Dipole 

For  a  slender  dipole  with  assumed  transverse  Maxwellian  distribution,  if  100  cells 
are  used  to  represent  the  current  longitudinal  dependence  and  50  points  are  used 
to  tabulate  the  function  5(A*),  then  in  order  to  fill  in  the  impedance  matrix 
the  proposed  algorithm  (algorithm  A)  requires  the  computation  of  an  infinite 
integral  150  times.  In  contrast,  the  algorithm  with  the  integral  in  polar  coordinate 
(algorithm  B),  typically  requires  the  computation  of  the  infinite  integral  more  than 
3000  times.  If  a  combination  of  the  piecewise  sinusoidal  and  Maxwellian  basis 
functions  are  used  in  a  Galerkin  procedure  of  the  moment  method,  with  algorithm 
A  it  takes  only  about  6  minutes  on  an  IBM  AT  computer  and  about  35  seconds 
on  a  Compaq-386  computer  to  generate  all  the  impedance  matrix  elements.  With 
algorithm  B  it  takes  about  one  hour  on  an  IBM  AT  and  about  6  minutes  and  30 
seconds  on  a  Compaq-386.  The  computer  codes  using  algorithm  B  for  testing  are 
based  on  the  one  discussed  in  [23],  which  is  about  2  times  more  efficient  than  that 
discussed  in  [20]. 


4.3.2  Microstrip  Open-end  and  Gap  Discontinuities 


For  the  microstrip  open-end  problem,  modeling  of  semi-infinite  lines  is  required. 
A  combination  of  semi-infinite  traveling  wave  modes  and  local  subdomain  modes 
'll  fruitful  and  can  be  modified  easily  to  apply  to  different  geometries  [21,23]. 
The  efficiency  of  the  present  algorithm  is  clearly  shown  in  these  open-end  and 
gap  computations,  where  the  mixed  type  of  basis  functions  are  used.  This  is 
because  if  only  one  transverse  basis  function  is  used  for  the  longitudinal  current, 
the  function  5 (A*)  is  the  same  for  all  the  impedance  matrix  elements.  Therefore, 
it  needs  to  be  computed  only  once.  The  computer  time  required  to  compute 
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the  equivalent  circuit  of  this  semi-infinite  microstrip  line  is  tested  for  both  the 
algorithm  with  a  polar  coordinate  integration  (algorithm  B)  and  the  one  with  a 
rectangular  coordinate  (algorithm  A)  integrsition.  For  a  Compaq-386  computer, 
with  20  subdomain  basis  functions,  algorithm  A  requires  about  30  seconds,  while 
algorithm  B  requires  about  200  seconds,  to  compute  one  data  point. 


4.4  Conclusions 

In  this  work,  a  new  algorithm  wzis  developed  for  the  method  of  moments  solution 
of  microstrip  antenna  and  circuit  problems.  The  algorithm  is  found  very  efficient 
and  can  be  applied  over  almost  the  entire  frequency  range.  The  efficiency  and  ver¬ 
satility  of  the  new  algorithm  are  illustrated  through  the  examples  of  a  microstrip 
dipole  and  a  microstrip  open-end  and  gap.  The  new  algorithm  can  also  handle  a 
variety  of  different  types  of  basis  functions,  which  will  find  many  applications  for 
the  analysis  of  irregular  shape  microstrip  antennas  and  circuit  discontinuities. 


Appendix  4-A 


An  example  is  presented  that  illustrates  the  integration  contour  deformation  that 
is  used  to  obtain  a  fast  convergent  integral.  Considering  two  x  roof-top  current 
basis  functions  reaction  with  the  pulse  functions  of  size  Wg  with  separation  y„  in 
a  microstrip  structure  with  substrate  thickness  h  and  dielectric  constant  £,,  for 
^  ^  and  A  »  ki,  5  (A*)  defined  in  Eq.  4.11  is  asymptotically  in  the  form  of 


S(A.)  =  P(A.)<3(A,) 

(4.12) 

where 

and 

P(A.)  = 

4  1  +  Cr 

(4.13) 

(4.14) 

Q{\x)  can  be  written  as 

<3(A.)  =  <?,(A.,<,)  -  i|Q,(A.,<ii)  +  (2,(A.,y 

(4.15) 

where 

(4.16) 

A 

II 

(4.17) 

1 

II 

(4.18) 

and 

^3  =  yr.  -1-  2u;e. 

(4.19) 

When  6  =  0,  Qi(A*,6)  has  an  exact  form.  When  6^0,  the  integration  contour 
can  be  deformed  in  the  complex  domain  along  the  contour  of  Re  (Ay  =  A)  and 
Im{Ay  >  0).  The  final  form  of  Qi(Ax,6)  to  be  evaluated  numerically  is 

e.(A..<)  =  | 

where 

/w - ,  ’ 

{A6  +  jxYyJ  (A,6)2  +  (A6  +  jxY 

(4.20) 

(4.21) 
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Chapter  5 


Infinite  Arrays  of  Microstrip 
Antennas  on  Generalized 
Anisotropic  Substrates 

In  this  chapter  the  problem  of  infinite  arrays  of  printed  antennas  on  generalized 
anisotropic  substrates  is  addressed.  The  antenna  elements  studied  include  mi¬ 
crostrip  dipoles  and  probe-fed  patches.  In  the  analysis,  the  substrate  anisotropy  is 
general  in  form,  with  nine  components  in  both  the  permittivity  and  permeability 
tensors.  This  allows  for  the  study  of  antennas  on  uniaxial  substrates  with  tilted 
optical  axis  or  on  gyrotropic  substrates.  A  rigorous  full-wave  moment  method  solu¬ 
tion  is  adopted.  Numerical  results  indicate  that  by  neglecting  substrate  anisotropy, 
the  performance  of  printed  antennas  may  be  greatly  affected.  Results  include  cases 
of  biased  ferrite  substrates,  in  which  the  antenna  characteristics  may  be  changed 
(or  controlled)  by  varying  the  bias  magnetic  fields. 


5.1  Introduction 

The  material  properties  of  the  substrate  significantly  affect  printed  antenna  perfor¬ 
mance  [1,2,27].  Many  practical  microwave  substrates  exhibit  anisotropy.  There  are 
two  types  of  material  anisotropy,  electric  and  magnetic  anisotropy.  For  materials 
with  electric  anisotropy,  the  permittivity  may  be  different  in  different  directions. 
Along  each  of  the  three  principal  axes  the  permittivity  may  be  characterized  by 
a  different  constant.  Useful  microwave  substrates  with  electric  anisotropy  include 
single-crystal  sapphire  (f„  =  =  9.4,  =  11.6),  boron  nitride  (f,,  =  = 

5.12,  =  3.4),  and  Epsilam-10  (t„  =  =  13,  e„  =  10.2).  Magnetized  ferrites 
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exhibit  magnetic  anisotropy,  where  the  permeability  is  in  the  form  of  a  third-rank 
tensor.  The  elements  of  the  permeability  tensor  are  related  to  the  external  applied 
dc  magnetic  field,  the  microwave  frequency,  as  well  as  the  inherent  physical  prop¬ 
erties  of  the  ferrite  material.  These  ferrite  materials,  which  are  usually  gyrotropic, 
are  tiseful  in  both  active  and  passive  circuits  due  to  the  fact  that  the  permeability 
tensor  and  thus  the  electrical  properties  of  the  material  are  changed  externally  by 
way  of  the  applied  dc  magnetic  field.  This  aspect  may  lead  to  many  applications 
in  printed  circuit  antennas  [28,29]. 

Infinite  array  analysis  is  a  useful  tool  for  the  design  of  large  antenna  arrays.  It 
provides  a  means  of  characterizing  some  of  the  important  features  of  the  large 
finite  array,  including  the  variation  of  element  input  impedance  as  a  function  of 
scan  angle.  The  use  of  dielectric  cover  layers  to  improve  the  input  impedance 
match  has  been  described  [32,33,34].  Infinite  array  analysis  has  been  applied  to 
microstrip  phased  arrays  on  isotropic  dielectric  substrates  by  several  investigators 
[30]-[34].  Pozar  and  Schaubert  discussed  the  phenomenon  of  scan  blindness  in 
infinite  phased  arrays  of  printed  dipoles  [30].  They  also  studied  the  infinite  array 
of  patches  [31].  The  work  reported  here  represents  a  significant  extension  of  those 
analyses  to  the  case  of  generalized  anisotropic  substrates. 

The  cited  work  employed  full-wave  method  of  moments  solutions  [30]-[34].  That 
formulation  has  been  shown  to  be  accurate  and  efficient.  A  fullwave  Green’s  func¬ 
tion  and  method  of  moments  integral  equation  solution  for  the  anisotropic  is  pre¬ 
sented  here.  Unlike  in  the  case  of  the  isotropic  substrate,  the  Green’s  function  for 
a  generalized  anisotropic  substrate  is  not  a  simple  closed  form.  The  exponential- 
matrix  method  [35]  provides  a  useful  technique  for  the  generation  of  the  spectral 
dyadic  Green’s  function. 

The  most  difficult  aspect  of  the  analysis  of  microstrip  antenna  problems  is  the 
modelling  of  the  feeds.  For  printed  dipoles,  a  delta-gap  generator  is  usually  as¬ 
sumed  to  simulate  the  two-wire  transmission  line  feed.  For  probe-fed  patches  the 
probes  can  be  modelled  as  ideal  current  generators.  This  has  been  found  to  be  a 
very  reasonable  model  for  thin  substrate  cases  [31]. 

In  this  chapter  the  effects  of  substrate  anisotropy  on  the  scan  characteristics  of 
infinite  phased  arrays  of  printed  dipoles  and  patches  are  investigated.  This  first 
centers  on  the  study  of  the  uniaxial  substrate  with  arbitrarily  oriented  optical  axis, 
in  which  case  all  nine  components  of  the  permittivity  tensor  are  nonzero.  Then 
the  focus  is  shifted  to  the  study  of  the  scan  characteristics  of  antennas  on  a  ferrite 
substrate.  The  effects  of  the  applied  dc  magnetic  field  and  its  direction  on  the 
array  performance  are  illustrated  through  numerical  examples. 

In  Section  5.2  the  infinite  array  analysis  including  both  permittivity  and  perme- 
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ability  anisotropy  is  described.  The  Green’s  function  for  a  generalized  anisotropic 
substrate  is  formulated,  where  both  horizontal  and  vertical  current  sources  are 
included. 

In  Section  5.3  the  importance  of  substrate  anisotropy  is  highlighted  through  a 
number  of  numerical  examples.  Some  features  of  infinite  arrays  of  printed  antennas 
on  a  ferrite  substrates  are  also  discussed. 


5.2  Analysis 


The  geometries  of  infinite  phased  arrays  of  printed  dipoles  and  patches  are  shown 
in  Figures  5.1  and  5.2,  respectively.  The  substrate  thickness  is  d.  All  the  dipoles 
are  identical  with  length  /  and  width  w.  All  the  patches  are  identical  with  length 
Lx  and  width  Ly.  The  structure  is  periodic  and  the  Floquet  cell  dimensions  are 
a  and  6  in  the  x-  and  y-  directions,  respectively.  The  antenna  elements  lie  in  the 
x-y  plane.  The  material  constants  of  the  grounded  substrate  are  assumed  to  have 
the  generalized  form: 


t^xy 

/ix* 

Cxx 

C*y 

Cx.  ' 

Mv* 

fiyy 

Mvx 

and  c  =  £o 

€yx 

£yy 

Cy* 

Mxx  . 

e.y 

(5.1) 


In  progressively  phased  arrays,  the  current  on  each  antenna  must  be  phased  in 
accordance  with  the  scan  angle  {6,  <!>)  as 


where 


=  sin^  cos^,  ’9/y  =  sin^  sin<^  , 


and  m  and  n  are  the 


indices  for  each  antenna  in  the  infinite  array. 


(5.3) 


The  Green’s  function  for  an  infinite  array  of  printed  antennas  [30]  has  the  form: 

G(^.  S')  =  A  E  E  <5(*-  *v)  e'*-'— 'I  .».(■'-*')  (5.4) 

m=:-oo  n=-cx> 


where 


kx 


2irm 

a 


+  ko'^z, 


ky 


2nn 

~r 


+  ko9y 


(5.5) 


and  G{kx,  ky)  is  the  spectral  Green’s  function.  In  the  case  of  the  generalized 
anisotropic  substrate,  this  spectral  Green’s  function  may  be  obtained  by  using  a 
exponential  matrix  method. 
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Figure  5.1:  Infinite  Array  of  Printed  Dipoles  on  a  Generalized  Anisotropic  Sub¬ 
strate. 


Figure  5.2:  Infinite  Array  of  Probe-fed  Patches  on  a  Generalized  Anisotropic  Sub¬ 
strate. 
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(5.6) 


Using  the  2-D  Fourier  transformation 

€  =  i  v  t  i 

H  ab  ^  ^  I 

n=-oo  m=-oo 


(ik.,  ky)  «-'***  e"^***'  dk^  dky, 


in  Maxwell’s  curl  equations,  one  obtains 


-  jk  X  ^  +  —{z  xH)  =  jui-  E  +  J6{z-  z') 
oz 


-  jk  X  .E  +  —  (2  X  E)  =  -jwA  •  H,  (5.8) 

az 

with  k  =  k^x+kyy.  E  and  H  are  the  spectral  electric  and  magnetic  fields;  moreover 

E  =  E,i  +  E„y  +  EgZ  (5.9) 

and 

H  =  H,x  +  Hyy  +  HJ.  (5.10) 

After  some  algebraic  manipulations,  Eqs.  5.7  -  5.8  can  be  written  in  the  form  of 
a  set  of  linear  differential  equations: 


—  |^(r)|  =  |.4||ii(z)]  +  |/|<(z-y) 


(5.11) 


where 


(^(2)]  = 


k,H,{z)  +  kyHy{z)  ■ 

kyHg{z)  —  kgHy{z)  ,  ,  ,  _ 

kgE,{z)  +  k„E„(2) 

A:„E*(2)  -  kgEy{z)  , 


an  an  an  au 

®21  ®22  ®23  ®24 

O31  O32  O33  <*34 

021  <142  <*43  <*44 


(5.12) 


In  Eq.  5.11,  [rp]  isa4xl  vector  containing  the  Fourier  transforms  of  the  tangential 
E  and  H  fields.  [A]  is  a  4  x  4  complex  matrix  whose  elements  are  completely 
described  by  the  material  parameters  f  and  and  the  Fourier  transform  variables 
k*  and  ky.  [/]  is  a  4  x  1  vector  containing  the  Fourier  transforms  of  any  sources 
that  might  be  in  the  layer.  For  horizontal  current  sources  one  has 


I/I  =  |t/.i  = 


or  i/I  =  |t/,l  = 


(5.13) 


[Ug\  is  for  the  current  source  in  the  x  direction,  while  \Uy\  is  for  the  current  source 
in  the  y  direction.  For  vertical  current  sources  (probes)  one  writes 


I/I  =  IC^.I  = 


0 


(5.14) 
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The  solution  of  Eq.  5.11  can  be  found  by  using  the  Cayley-Hamilton  theorem  [35] 
or  using  eigenvector  analysis.  The  results  from  either  method  are  equivalent  and 
yield  terms  each  of  which  has  a  z-dependence  of  the  form  where  A,  are  the 
eigenvalues  of  the  [A]  matrix. 


5.2.1  Integral  Equation  and  Green's  Function  for  Infinite 
Probe-fed  Patch  Arrays 


The  integral  equation  for  infinite  probe-fed  patch  arrays  has  the  form 


Gx*  G,, 

GyX  Gyy 


dsi  -f 


Jtds2  =  0, 


(5.15) 


where  G^y . are  in  the  spectral  form  of  Eq.  5.6. 


To  find  the  dyadic  Green’s  function  for  patch  antennas,  we  describe  the  fields 
by  a  vector  [tp]  given  by  Eq.  5.12.  The  unknown  constants  can  then  be  found  by 
applying  the  appropriate  boundary  conditions  at  each  layer  interface.  The  solution 
of  Eq.  5.11  is 

[0(2i)j  =  [T{zi  -  Z2)]  [^{Zi)],  (5.16) 


where 


■  e*'*  0  0  0' 

|r(^)l  =  W  °  'o'  X-  0  I*!"' 

0  0  0 


(5.17) 


in  which  the  A,,  t  =  1,2, 3,4  are  the  eigenvalues  and  [<^j  is  a  4  x  4  eigenvector 
matrix  of  [Aj.  (For  the  isotropic  case  there  would  be  repeated  eigenvalues,  and 
the  terms  of  the  exponential  matrix  could  be  greatly  simplified  and  would  exhibit 
sink  and  cosh  behavior.)  According  to  the  Calley-Hamilton  theorem,  the  matrix 
[r(2)]  can  also  be  expressed  as 


|r(z)|  =  E 


(5.18) 


where  the  coefficients  C„  are  taken  from  the  solution  of  the  matrix  equation 


.W  =  ^A‘C»{z),  i=l,2,3,4. 


(5.19) 


It  is  quite  clear  that  the  physical  meaning  of  the  matrix  [Tj  is  to  relate  the  tan¬ 
gential  electromagnetic  fields  at  one  surface  {z  =  Zi)  to  another  {z  =  Z2).  For 
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the  geometry  shown  in  Figure  5.2,  the  tangential  electromagnetic  fields  at  the 
air-ferrite  interface  can  now  be  related  to  those  at  the  ground  plane: 


[^(d)]  =  [T{d)]  [0(0)] . 


(5.20) 


The  electromagnetic  fields  in  the  air  {z  >  d)  can  be  derived  in  a  straightforward 
manner.  With  this  result  together  with  the  boundary  conditions  that  the  tangen¬ 
tial  electric  fields  are  zero  at  the  ground  plane,  one  obtains 


(5.21) 


jy/k’  -  kl  a  '  c  ■ 

1  "  n  ’ 

3\Jk?-klb  u 

— (j/xo  a  .  L  ®  . 

where  d,  6,  c  and  d  are  unknown  spectral  quantities  to  be  determined. 


With  the  current  source  at  the  air-ferrite  interface,  the  electromagnetic  fields  just 
above  and  below  z=d  are  related  as 


[0(d)]  -  [r(d)][0(o)]  =  [UjI 


(5.22) 


where  [Ui]  is  either  [l/*|  or  [Uy].  Eq.  5.22  involves  a  4  by  4  matrix  equation,  where 
the  solution  determines  the  quantities  d,  6,  c  and  d.  As  a  result,  the  spectral 
Green’s  functions  for  patches  are  determined  from  the  expressions 


G’ii(/Ci,  fcy)  or  Gxy(fcxi  kyj  — 


jkxJkl  +  A:*  -  k^b  -  unokyd 


kl  +  kl 


(5.23) 


^yxfykx,  ky)  or  Gyy{^kx,  ky^  — 


jkyyjkl  +  kl  -  k^b  +  ujfiokxa 


(5.24) 


When  |17/]  =  \Ux],  the  above  two  equations  are  for  G!xx{kx,  ky)  and  Gyx{kx,  ky); 
while  when  [Uj]  =  \Uy],  the  above  two  equations  are  for  Gxy{kx,  ky)  and  Gyy{kx,  ky). 

The  Green’s  function  for  the  probes  can  be  found  by  solving  the  following  matrix 
equation. 

\T{z‘  -  rf)||^(rf)l  -  |r(r’)||lj(0)|  =  lU.I,  (5.25) 


|0(i)|  -  (r(<i)|10(O)l  =  {T{d  -  z')l|f/.l,  (5-26) 

where  [17,]  is  defined  in  Eq.  5.14.  With  [0(d)]  and  [0(0)]  defined  in  Eq.  5.21,  the 
expression  for  G,,  and  Gy,  can  be  found  from  Eqs.  5.23  and  5.24  by  replacing  G,, 
by  G„  and  G,y  by  Gy,. 
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5.2.2  Method  of  Moments  for  Infinite  Patch  Arrays 


In  the  moment  method  procedure,  the  unknown  current  distribution  on  the  patches 
is  expanded  in  terms  of  a  set  of  entire  domain  basis  functions: 


J.  =  '£rt  y) 

i 

(5.27) 

and 

!/). 

• 

(5.28) 

where 

Jf  =  ^  sin  ^(x -h  L,/2)  cos  ^(y -h  Ly/2) 

X/y 

(5.29) 

and 

1  ir'tr  Itt 

Ji  =  —  cos  —  (x  +  i/,/2)sin  —  (y -t- Ly/2). 

Xy*  Ijti 

(5.30) 

The  current  on  the  probes  is  assvuned  to  be  constant  along  the  vertical  direction 
(z  axis),  and  to  have  a  delta-function  form  in  the  x-y  plane:  6{x  —  Xp)6{y  —  j/p). 


With  the  same  expansion  and  testing  functions  in  the  integral  equation,  one  ob¬ 
tains  the  matrix  equation 

X;  Zijij  =  Vi  for  all »,  (5.31) 

i 

where  is  in  the  form  of 

2'.  =  i  E  E  •'<'5  j,*  .  (5.32) 

m  n 

G  being  any  of  the  four  dyadic  components  G**,  G*y,  Gy,,  Gyy,  and  J,  ,  the  Fourier 
transform  of  the  expansion  or  testing  function.  The  components  of  the  excitation 
or  voltage  matrix  which  is  due  to  the  probe  current  take  the  form  of 

=  ^EEj/G„e--‘-'^e-'^  (5 

m  n 

or 

=  4E  (5.34) 

m  n 

The  method  of  moments  provides  the  currents  on  the  patches.  The  input  impedance 
from  the  probe  is  calculated  as 

^  E^.dz,  (5.35) 

I*  Jo 
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where  Ip  is  the  current  on  the  probe  and  E^{z)  is  the  electric  field  due  to  the 
currents  on  the  patch.  If  reciprocity  applied,  the  input  impedance  could  be  written 

=  -HW’  (5-36) 

3 

where  Vj  is  but  the  matrix  element  of  the  excitation  matrix  in  Eq.  5.33.  However, 
for  the  generalized  anisotropic  substrate  reciprocity  no  longer  holds.  Alternative 
procedures  for  the  computation  of  the  electric  fields  due  to  the  patch  currents  are 
required.  One  can,  from  Eq.  5.22,  find  [^(0)]  and  obtain  the  lateral  electromag¬ 
netic  fields  at  any  position  z  from  the  expression 

|v;{^)l  =  |r(^)|  |t^(0)l.  (5.37) 

Further,  the  spectral  electric  field  can  be  related  to  the  lateral  electromagnetic 
fields  through 

f^oitzE^  =  kyHj,  -  ktHy  -  f-ottzEi  -  eoe^yEy.  (5.38) 

Replacing  and  G^y  in  Eq.  5.33  by  this  spectral  field  E,,  one  obtains  a  voltage 
matrix  Vj’.  As  a  result,  the  input  impedance  of  the  probe  is  calculated  as 

^in  =  -EW-  (5-39) 


5.2.3  Integral  Equation  and  Green’s  Function  for  Infinite 
dipole  Arrays 


For  narrow  microstrip  dipoles,  the  transverse  current  is  usually  negligible.  There¬ 
fore,  the  corresponding  integral  equation  is  simplified  to  the  form 


J  a zx  'It  3?a)) 


(5.40) 


where  Xa  is  the  location  of  the  delta  gap  source.  The  spectral  form  of  Gn  is  the 
same  as  the  one  discussed  previously  for  patches.  If  roof-top  basis  functions  are 
used  in  the  method  of  moments,  the  impedance  matrix  elements  can  be  expressed 
ais 


ab 


f;  f;  G{k.,  k.)  j\k.,  *,) 

m=-oo  n=-oo 


(5.41) 


where 


.  sin(fcyty/2)  1  —  cos(fciA) 

>  - —  , 


kyW 


(5.42) 


and  A  is  half  the  length  of  the  basis  function.  For  a  center-fed  dipole,  it  is 
convenient  to  use  an  odd  number  of  basis  functions  such  that  all  the  elements 
in  the  voltage  vector  are  zero  except  the  center  term  which  is  unity.  The  input 
impedance,  which  is  the  reciprocal  of  the  current  at  the  center  of  the  dipole,  can 
be  obtained  from  the  matrix  equation  \Z]\I]  =  [V']. 
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5.3  Results 


The  results  shown  first  illustrate  the  effects  of  substrate  anisotropy  on  the  scan 
chau'acteristics  of  infinite  phased  arrays  of  printed  dipoles.  For  a  uniaxial  substrate 
with  an  optical  axis  tilted  in  the  x-z  plane,  where  the  optical  axis  angle  measured 
from  the  z  axis  is  Bq,  the  permittivity  tensor  is 


6  =  to 


e,  cos^  Bq  -f  f,  sin’  Bq  0  (c*  —  e,)  sin  Bq  cos  Bq 

0  c.  0 

I  (ti  —  t,)  sin  Bq  cos  Bq  0  €*  sin’  +  t,  cos’  Bq 


(5.43) 


Figures  5.3  and  5.4  show,  respectively.  E-plane  and  H-plane  scan  characteristics, 
comparing  cases  with  and  without  material  anisotropy.  The  substrate  is  PTFE 
cloth,  which  exhibits  uniaxial  anisotropy  [36j.  It  is  seen  that  the  effects  of  substrate 
anisotropy  are  quite  noticeable  for  the  E-plane  but  not  for  the  H-plane.  Also,  the 
anisotropy  in  the  planar  direction  (x-y  plane)  is  more  important  than  that  in  the 
vertical  direction.  Figures  5.5  and  5.6  show  the  scan  characteristics  for  the  E  and 
H-planes,  respectively,  for  a  sapphire  substrate.  The  characteristic  impedance  of 
the  feed  is  assumed  matched  with  the  input  impedance  of  the  dipole  at  a  broadside 
scan,  and  when  the  substrate  is  isotropic  with  e,  =  9.4.  In  addition  to  the  results 
for  the  isotropic  case,  the  results  for  four  different  directions  of  the  optical  axis 
are  also  shown  in  Figures  5.5  and  5.6  (5o  =  0°,  30®,  60°,  90°).  It  is  seen  that  the 
substrate  anisotropy  introduces  a  severe  input  mismatch,  although  the  angle  of 
scan  blindness  does  not  change  much.  The  results  also  show  that  when  Bq  =  0° 
the  mismatch  is  most  severe.  In  other  words,  the  anisotropy  in  the  z  direction  is 
more  important  than  that  in  the  x  direction. 


Examples  of  scan  characteristics  for  an  infinite  array  of  probe-fed  patches  on 
Epsilam-10  substrate  are  shown  in  Figures  5.7  and  5.8  for  the  E  and  H-planes, 
respectively.  The  results  for  the  optical  axis  in  Bq  =  0°,  30°,  60°,  90°  are  shown. 
It  is  seen  again  that  the  material  anisotropy  introduces  large  changes  in  the  patch 
array  element  input  impedance.  Ignoring  the  material  anisotropy  in  such  a  case 
would  result  in  a  severely  mismatched  design. 


The  permeability  tensor  of  a  ferrite  substrate  is 


M  = 


H+  {1  -  n)  cos^  4>o 
{1  —  fi)  sin(^ocos<^o 
-j/csin  <j>o 


(1  —  A*)  sin(^ocos<^o 
Ai  +  (1  -  A*)  sin’(;io 
Jk  cos  <f>o 


j/csin  <f>o 
-jK  cos  <po 


(5.44) 


where  4>q  is  the  angle  of  the  direction  of  the  dc  bias  H  field  measured  from  the  x  axis 
in  the  x-y  plane.  The  parameters  //  and  k  are  characterized  by  the  dc  magnetic 
field,  the  saturation  magnetization,  and  the  frequency  of  operation  [9].  The  effect 
of  the  applied  magnetic  field  on  infinite  printed  dipole  arrays  is  shown  in  Figure 
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Figure  5.3:  Comparison  of  E-plane  Scan  Characteristics  of  Infinite  Printed  Dipole 
Array  on  Uniaxial  and  Isotropic  Dielectric  Substrates. 

Reflection  coefficient  magnitude  for  the  E-plane 

d  =  0.19A„,  a  =  0.5155A„,  b  =  0.5A„,  L  =  0.39Ao,  w  =  0.002A„. 
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Figure  5.4:  Comparison  of  H-plane  Scan  Characteristics  of  Infinite  Printed  Dipole 
Array  on  Uniaxial  and  Isotropic  Dielectric  Substrates. 

Reflection  coefficient  magnitude  for  the  H-plane  scan. 


d  =  0.19A,,  a  =  0.5155Ao,  b  =  0.5A„,  L  =  0.39A„,  w  =  0.002Ao. 
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Figure  5.5:  E-plane  Scan  of  a  Dipole  Array  on  a  Uniaxial  Substrate  with  Tilted 
Optical  Axis 

Reflection  coefRcient  magnitude  for  the  E-plane  scan. 

The  optical  axis  tilt  angle  is  Og. 

Zin  =  8.06  +  jO.36. 

£,  =  £„  =  9.4,  €g  =  11.6. 

d  =  0.06Ap,  a  =  0.5Ao,  b  =  0.5Ao,  L  =  O.I86A0,  it;  =  0.002A„. 
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Scan  Angle  in  Degree 

Figure  5.6:  H-plane  Scan  of  a  Dipole  Array  on  a  Uniaxial  Substrate  with  Tilted 
Optical  Axis 

Reflection  coefficient  magnitude  for  the  H-plane  scan. 

The  optical  axis  tilt  angle  is  Og. 

Zin  =  8.06  +  jO.36. 
c,  =  e„  =  9.4,  e,  =  11.6. 

d  =  0.06A„,  a  =  0.5Ao,  6  =  0.5A<„  L  =  0.186A<„  w  =  0.002Ao. 
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con  angle  in  degrees 


Figure  5.7:  E-plane  Scan  Characteristics  of  an  Infinite  Probe-fed  Patch  Array  on 
a  Uniaxial  Substrate  with  Tilted  Optical  Axis. 

Reflection  coefficient  magnitude  for  the  E-plane  scan. 

The  optical  axis  tilt  angle  is  Bo.  f  =  4.085GHz. 

Lg  =  1.1cm,  Ly  =  1.1cm,  Xp  =  0.165cm,  j/p  =  0.0cm, 
e,  =  =  13.0,  £,  =  10.2. 

d  =  0.127cm,  a  =  2.125cm,  6  =  2.378cm,  =  47  fl. 
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Figure  5.8:  H-plane  Scan  Characteristics  of  an  Infinite  Probe-fed  Patch  Array  on 
a  Uniaxial  Substrate  with  Tilted  Optical  Axis. 


Reflection  coefficient  magnitude  for  the  H-plane  scan. 
The  optical  axis  tilt  angle  is  /  =  4.085GHz. 

Lx  =  1.1cm,  Ly  =  1.1cm,  Xp  =  0.165cm,  yp  =  0.0cm, 

C,  =  =  13.0,  £x  =  10.2. 

d  =  0.127cm,  a  =  2.125cm,  6  =  2.378cm,  Zc  =  47  Cl. 
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Figure  5.9:  Comparison  of  the  Scan  Characteristics  of  Dipole  Arrays  on  Isotropic 
and  Biased  Ferrite  (gyrotropic)  Substrates. 

d  =  O.O6A0,  a  =  O.5A0,  6  =  O.5A0,  Cr  =  12.6,  u;  =  0.002Ao. 

For  the  isotropic  substrate  L  =  0.156Ao  and  Zc  =  9.3  +  j0.30. 

For  the  ferrite  substrate  L  =  0.205Ao  and  Zc  =  0.7  +  j0.03. 


Figure  5.10:  E-plane  Scan  Characteristics  of  an  Infinite  Array  of  Printed  Dipoles 
on  a  Ferrite  Substrate. 

d  =  O.OeA,,,  a  =  0.5Ao,  b  =  O.SA^,  =  12.6,  w  =  0.002Ao. 

The  dc  magnetic  field  direction  is  M  = 


4>q  =  0°,  I  =  0.205Ao,  Zc  =  0.7  +  j0.03 
<^o  =  30°,  L  =  0.1893Ao,  Zc  =  0.96  -  j0.02 
<i>o  =  90°,  L  =  0.1603Ao,  Zc  =  10.2  -  yo.06 


Scan  Angle  in  Degree 


Figure  5.11:  H-plane  Scan  Characteristics  of  an  Infinite  Array  of  Printed  Dipoles 
on  a  Ferrite  Substrate. 


d  =  O.O6A0,  a  =  O.SAo,  b  =  0.5Ao,  c,  =  12.6,  u;  =  0.002Ap. 

The  dc  magnetic  field  direction  is  (^<,,^0),  M  =  0.7665,  k  =  —0.5665. 
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Figure  5.12:  Scan  Characteristics  of  an  Infinite  Probe-fed  Patch  Array  on  a  Ferrite 
Substrate. 

/  =  0.66GHz,  d  =  0.3175cm, 

Lz  =  2.98cm,  Ly  =  2cm,  Xp  =  0.88cm,  yp  =  0.0cm, 

Cr  =  12.6,  HoHo  =  0.05T,  MoM,  =  0.275T,  6^  =  90",  4>o  =  0", 
a  =  5  cm,  b  =  4  cm,  and  Ze  =  50Q. 
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Figure  5.13:  Input  Reflection  Coefficient  for  the  Patch  Array  on  a  Biased  Ferrite 
Substrate  as  a  Function  of  the  Biais  Field  Angle. 

/  =  0.66GHz,  d  =  0.3175cm, 

Lx  =  2.98cm,  L„  =  2cm,  Xp  =  0.88cm,  =  0.0cm, 

er  =  12.6,  ^loHo  =  0.05T,  ^loM.  =  0.275T,  =  90°, 

a  =  5  cm,  b  =  4  cm,  and  Ze  =  SOfl. 
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5.9.  The  results  for  an  isotropic  substrate  with  no  bias  are  also  shown.  The  results 
for  a  ferrite  substrate  are  for  the  case  in  which  the  applied  magnetic  field  is  in  the  x 
direction  (same  as  the  current).  It  is  seen  that  the  input  impedance  level  reduces 
drastically  with  a  ferrite  substrate.  The  results  also  show  that  the  scan  angle 
where  the  scan  blindness  occurs  can  be  much  larger  (>  80“)  with  a  biased  ferrite 
substrate  (achieved  without  a  reduction  of  the  Floquet  cell  size).  The  effects  of  the 
direction  of  the  bias  magnetic  field  on  the  scan  characteristics  are  shown  in  Figures 
5.10  and  5.11  for  the  E  and  H-planes,  respectively.  It  is  seen  that  when  the  dc 
magnetic  field  direction  is  transverse  to  the  direction  of  the  current  in  the  dipole, 
the  scan  characteristics  are  closer  to  the  isotropic  case  (no  applied  dc  magnetic 
field).  When  the  direction  of  the  dc  magnetic  field  decreases  from  =  90°,  the  E- 
plane  performance  improves,  while  the  H-plane  scan  characteristics  are  degraded. 

The  scan  characteristics  for  an  infinite  probe-fed  patch  array  on  a  biased  ferrite 
substrate  are  shown  in  Figure  5.12.  It  is  seen,  as  for  the  dipole  case,  that  the 
E-plane  and  H-plane  scan  characteristics  are  very  much  the  same.  However,  it  is 
found,  for  the  patch  array,  that  a  scan  blindness  condition  arises  in  the  H-plane 
but  not  in  E-plane.  Also  it  is  seen  that  the  scan  performance  is  much  better  in 
the  diagonal  plane  (^  =  45°),  where  even  at  the  extreme  scan  angle  of  80  degrees 
the  mismatch  not  too  large  ([Rl  <  0.40). 

It  is  known  that  the  bias  magnetic  field  changes  the  permeability  and  thus  the 
electrical  properties  of  ferrites.  It  is  therefore  possible  to  adjust  the  bias  H  field 
level  and  thereby  effect  changes  in  the  scan  performance  of  an  array.  It  is  also 
seen  from  Eq.  5.44  that  the  permeability  tensor  of  ferrites  is  a  strong  function  of 
the  angle  of  the  bias  H  field.  An  example  of  the  effect  of  the  bias  angle  on  the 
broadside  scan  input  impedance  is  shown  in  Figure  5.13.  As  <f>o  varies  from  0  to 
90  degrees,  the  input  impedance  changes  from  that  of  a  perfect  match  to  one  of 
almost  total  reflection. 


5.3.1  Conclusion 

An  analysis  of  infinite  phased  arrays  of  printed  antennas  on  generalized  anisotropic 
substrates  was  developed.  It  has  been  applied  to  uniaxial  dielectric  substrates,  and 
to  gyrotropic  (biased  ferrite)  substrates.  From  the  computed  results,  it  was  found 
that  the  substrate  dielectric  anisotropy  in  commonly  used  substrates  can  introduce 
significant  changes  in  active  impedance  of  the  array  element.  If  those  anisotropic 
effects  are  neglected  in  the  design  process  the  performance  of  the  array  can  be 
severely  degraded.  In  the  case  of  a  ferrite  substrate  it  was  also  found  that  the 
strength  and  the  direction  of  the  applied  magnetic  field  may  affect  significantly 
the  scan  characteristics  of  the  infinite  array.  It  was  observed,  in  the  principal 
planes,  that  the  angle  at  which  scan  blindness  occurs  may  be  very  large  when  a 
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biased  ferrite  substrate  is  used,  and  that  this  condition  can  be  aichieved  without 
reducing  the  size  of  the  periodic  cell. 
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